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ABSTRACT 

;  A  UHF  clear-air  Doppler  radar  wind  profiler  was  purchased  by  the  Naval  PoSt- 

I  graduate  School  and  began  operation  in  late-March  of  1990.  The  identification  of  nu- 

t 

merous  meteorological  phenomena  utilizing  UHF  Doppler  radar  wind  profilers  has  been 
1  successfully  documented;  however,  published  work  to-date  has  concentrated  on  inter- 

I  pretation  of  hourly-consensus  data.  But  with  the  onset  of  improved  technology,  it  is 

I  «  now  possible  to  retrieve  reliable  data  within  the  hourly  period.  This  study  examines  the 

quality  of  the  six-minute  data  through  statistical  analysis  and  determines  the  overall  ac- 
-  curacy  of  the  School's  profiler.  Specifically,  the  performance  of  the  profiler  is  compared 

with  rawinsonde  data,  other  profiler  comparison  studies,  as  well  as  the  contract  specifi¬ 
cations.  'The  study  also  looks  at  overall  performance  of  UHF  Doppler  radar  profiler  in 
an  environment  characterized  by  a  strong  marine  inversion.  Finally  the  study  presents 
evidence  that  meteorological  phenomena  which  exhibit  time  scales  on  the  order  of  min¬ 
utes  can  be  observed  using  the  six-minute  profiler  data. 
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I.  INTRODUCTION 


Radar  technology  has  undergone  continuous  refinement  since  it  was  first  developed 
early  in  this  century.  The,  development  of  the  modem  wind  profiler  is  an  outgrowth  of 
research  done  in  the  I920's  and  I930's  with  radars  designed  to  probe  the  ionosphere. 
This  research  was  driven  by  the  need  to  calculate  HF  radio  interference,  for  example, 
from  the  velocity  of  the  free-electrons  and  ions.  Interest  faded  through  the  late  1930's. 
It  was  not  until  World  War  II  that  work  began  again  with  short-wave  radars-hoping 
again  for  the  chance  to  utilize  the  ionosphere  successfully  for  communications.  Al¬ 
though  the  initial  applications  of  radar  were  for  communications  research  and  detecting 
solid  objects,  such  as  warships,  it  was  soon  discovered  that  radar  could  also  detect 
hydrometeors  in  the  atmosphere.  It  was  not  until  the  early  1970's  that  strong  returns 
were  first  received  in  "clear  air"  (Woodman  and  Guillen,  1974), 

In  the  late  1970's  the  first  significant  work  with  profilers  began  in  the  United  States. 
While  the  early  wind  profilers  were  designed  to  sample  the  stratosphere  and  even  the 
mesosphere,  today's  versions  concentrate  on  profiling  the  troposphere  and  lower 
stratosphere.  The  National  Oceanic  and  Atmospheric  Administration  (NOAA)  Envi¬ 
ronmental  Research  Laboratory  in  Boulder,  Colorado  began  research  with  50  MHz  and 
400  MHz  radars  in  the  early  1980's  and  since  then  have  been  operating  a  small  network 
of  wnd  profilers  in  Colorado  (see  inset  within  Fig.  1). 

For  the  past  seven  years  the  Denver  Weather  Service  Forecast  Office  has  had  access 
to  the  profiler  network  and  found  the  profiler  data  to  be  very  useful  in  forecasting  winter 
storms.  Although  the  most  experience  dealing  with  wind  profilers  and  their  subsequent 
network  capabilities  has  been  in  the  Colorado  area,  it  is  expected  that  profilers  will  be 
extremely  useful  in  other  climatic  regions. 

The  NOAA  Wave  Propagation  Laboratory  is  in  the  process  of  setting  up  a  demon¬ 
stration  network  of  30  wind  profilers  by  1991  (locations  depicted  in  Fig.  1  from  van  de 
Kamp  (1988)).  The  first  UNISYS-produced  profiler  was  officially  accepted  by  NOAA 
in  August  of  1989,  located  at  Platteville,  Colorado.  As  the  network  becomes  complete, 
forecasters  and  researchers  will  have  the  opportunity  to  compute  and  evaluate  hourly 
time  series  of  wind  profiles  and  derived  kinematic  quantities. 

The  Naval  Postgraduate  School  has  purchased  a  404  MHz  profiler  as  part  of  the 
Marine  Atmospheric  Measurements  Lab  in  Spanagel  Hall.  The  transceiver  antenna  was 
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Location  of  wind  profiltr  demonstrotion  network  sites 
to  be  operotionot  by  1991. 


Hg.  1.  Profiler  Sites  for  tlie  NWS  Demonstration  Network:  (Source:  van  dc 
Kamp{1988)), 

installed  at  Fritzclm  Field.  Fort  Ord.  The  area  topography  and  geographical  position 
ofthe  NPS  profiler  are  presented  in  Fig.  2.  Figure  3  portrays  details  ofthe  NFS  profiler 
site  and  flat-bed  antenna  array.  The  profiler  was  operational  by  late-March  1990. 

The  NPS  profiler  was  tlie  first  UMF  profiler  placed  on  the  California  central  coast, 
and  among  the  first  on  the  North  American  west  coast.  The  California  central  coast 
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Kilometers  east  of  Profiler  Site 


Fig.  2.  Topographical  Map  of  area  Surrounding  NFS  Profiler  Site 

provides  a  unique  environment  influenced  by  cold-water  upwelling,  a  persistent 
subsidence  of  dry  air  due  to  the  eastern  Pacific  subtropical  ridge,  and  the  complex 
interaction  of  topography  and  the  land/sea  breeze  circulation.  Consequently,  there  is 
need  to  determine  the  feasibility  of  operating  a  404  Mllz  ladar  profiler  in  an  environ¬ 
ment  characterized  by  a  strong  marine  inversion. 

Significant  work  with  profiler  data  has  already  been  accomplished,  especially  by  Kuo 
el  at.  (1987)  utilizing  the  Colorado  network  and  Forbes  ei  aL  (1989  and  1990)  with  the 
Penn  State  profilers.  Kuo  ei  al.  found  that  retrieved  temperature  derived  from  a  network 
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Fig.  3.  Configuration  of  the  NFS  Profiler  Site:  (Source:  Tycho  Technology,  Inc., 
Wind  Profiling,  1988) 


of  VHF  wind  profilers  was  significantly  more  accurate  than  the  direct  temperature 
measurement  from  a  combined  satellite-  and  ground-based  microwave  radiometric  sys¬ 
tem.  Forbes  et  al.  documented  how  wind  profiler  data  together  with  conventional 
National  Weather  Service  rawinsonde  data  proved  quite  helpful  in  understanding  the 
cause  of  subtle  weather  systems  seen  in  satellite  imagery.  Forbes  et  al.  also  examined 
the  utility  of  profiler-derived  temperature  gradients,  temperature  advections  and  lapse 
rate  gradient  information.  Both  the  Colorado  and  Pennsylvania  network  of  profilers 
have  been  operating  for  nearly  a  decade  and,  subsequently,  numerous  works,  beyond 
which  can  be  cited  in  this  study,  have  analyzed  important  aspects  of  their  respective  local 
meteorology,  including  frontal  passages  and  convective  storms. 

Published  studies  to-date  have  concentrated  on  interpretation  of  hourly-consensus 
data.  Forbes  et  al.  (1990)  discusses,  although  only  briefly,  interhour  fluctuations  de¬ 
tected  with  the  Penn  State  VHF  profiler,  but  their  emphasis  is  on  analyzing  significant 
mesoscale  weather.  The  primary  objective  of  this  study,  in  contrast,  is  examination  of 
the  potential  for  a  clear-air  wnd  profiling  Doppler  radar  to  observe  small-scale,  high- 
frequency  meteorological  cvcnls-in  other  words,  information  within  the  time  scale  of 
one  hour,  down  to  time  scales  of  six  nunutes  in  the  case  of  the  NPS  radar.  To  meet  this 
objective  it  will  be  necessarj’  to  first  ensure  data  from  the  individual  profiles  are  as  reli¬ 
able  as  possible,  i.e.,  throw  out  the  "bad"  data.  This  data  quality  process  will  hopefully 
provide  an  improved  capability  to  obtain  good  representative  data  beyond  that  which 
the  manufacturer's  software  provides.  It  is  also  hoped  that  by  improving  the  individual 
profiles  that  the  hourly-consensus  will  be  improved. 

A  secondary  objective  of  this  study  is  to  document  the  performance,  including  ac¬ 
curacy,  of  the  NPS  profiler  operating  in  an  environment  characterized  by  a  strong  ma¬ 
rine  inversion.  Determination  of  system  accuracy  is  important  in  this  study  in  order  to 
lend  insight  into  the  thermodynamic  relationships  presented  later.  The  study  will  also 
compare  the  data  accuracy  to  other  profiler  and  rawinsonde  studies. 

Finally,  this  study  hopefully  will  expose  many  students  and  instructors  to  the  capa¬ 
bility  of  a  UHF  radar  profiler,  in  addition  to  providing  future  students  the  necessary 
groundwork  to  expand  and  analyze  specific  meteorological  phenomena  with  high- 
temporal  resolution  and  on  a  real-time  basis. 

This  study  documents  research  performed  during  the  period  of  initial  operation,  30 
March  to  05  May  1990.  The  period  was  selected  for  several  reasons.  First,  the  profiler 
was  scheduled  to  be  part  of  Tropical  Cyclone  Motion  experiment  (TCM-90)  in  southeast 
Asia,  which  meant  disassembly  by  no  later  than  June  1990.  Second,  a  transmitter  am- 


5 


plifier  failed  one  month  prior  to  scheduled  disassembly,  and  the  amplifier  subsequently 
had  to  be  removed  and  sent  back  to  the  manufacturer.  The  period  left  (a  little  more 
than  a  month)  was  further  reduced  in  length  by  spurious  software  failures  and  power 
outages. 

Several  other  factors  limited  the  ability  to  conduct  a  detailed  meteorological  inves¬ 
tigation.  The  initial  operating  period  was  never  meant  to  be  analyzed  in  such  detail  as 
this  study  is  directed.  Therefore,  the  ancillary  data-gathering  systems  such  as  SODAR 
and  high-frequency  pressure  and  humidity  recording  devices,  were  never  arranged. 
Rawinsondes  were  the  only  equipment  configured  for  this  experiment,  but  at  the  last 
moment  the  primary  rawinsonde  system  allotted  for  comparison  studies  failed.  A  sec¬ 
ondary  unit  was  available  but  could  be  used  only  intermittently,  as  it  had  to  be  scheduled 
around  cruise  course  requirements.  In  the  end,  only  nine  successful  rawinsonde 
launchings  were  available,  and  five  of  those  were  from  onboard  the  student  cruise  boat 
■'•st  outside  Monterey  Bay. 

The  equipment  configuration  used  for  this  study  is  described  in  Appendix  A.  The 
Radio-Acoustic  Sounding  System  (RASS)  modification  for  the  wind  profiler  was  not 
installed  for  this  period,  and  the  software  for  Q  computations  was  not  yet  available 
from  the  manufacturer. 
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II.  BACKGROUND 


A.  THEORY  OF  PROFILER  OPERATION 

The  wind  profiler,  like  all  radars,  detects  the  transmission  or  reflection  of 
electromagnetic  radiation  at  certain  frequencies,  or  wavelengths.  A  phenomenon  known 
as  Bragg  scatter  is  critical  to  receiving  the  necessary  backscattered  signal  strength  in 
clear  air.  Bragg  scatter  comes  about  in  two  ways:  (I)  through  the  scatter  from  a  surface 
or  layer  with  a  wave  speed  which  is  resonant  to  the  cosine  of  the  incident  angle  of  the 
transmitted  radar  pulse,  or  (2)  through  the  presence  of  multiple  refractive  layers  spaced 
at  the  scales  of  turbulent  energy  as  is  in  the  case  of  UHF  profilerS“404  MHz  falls  in  the 
Ultra  High  Frequency  (UHF)  range.  Although  all  scales  of  turbulence  scatter  energy, 
it  is  not  until  the  turbulence  scales  approach  1/2  the  transmitted  wavelength,  or  37  cm 
in  the  case  of  the  404  MHz  profilers,  that  the  scattered  energy  waves  interact  construc¬ 
tively. 

Profilers  depend  upon  the  scattering  of  electromagnetic  energy  by  minor  irregulari¬ 
ties  in  the  index  of  refraction  of  air.  Energy  is  first  radiated  at  404  MHz  from  a  coaxial 
colincar  antenna  (Fig.  4),  essentially  several  cylindrical  patterns  radiating  from  inter¬ 
leaved  coaxial  cabling.  Above  approximately  200  m,  beyond  the  near-field,  these  indi¬ 
vidual  cylindrical  patterns  begin  to  interact  constructively  to  form  the  radio  wave  front. 
The  antenna  is  constructed  so  most  of  the  energy  is  embodied  within  a  five  degree  cir¬ 
cular  beam,  or  cone;  although  several  side  lobt  arc  present  as  illustrated  by  the  insert 
within  Fig.  4.  This  wave  front  is  not  impeded  intil  turbulence  on  the  correct  scale  is 
encountered,  as  shown  schematically  in  Fig.  4.  Turbulence  produces  strong  small-scale 
variability  in  three  dimensions,  carrying  moisture  and  temperature  discontinuities  verti¬ 
cally  (and  horizontally)  adjacent  to  other  air  parcels.  Atmospheric  turbulence  is  gener¬ 
ated  at  scales  of  tens  to  hundreds  of  meters  by  wnd  shear  and  convection.  These  large 
scale  eddies  break  up  into  smaller  eddies,  which  break  up  into  still  smaller  ones,  and  so 
on  until  the  turbulent  energy  is  dissipated  at  very  small  scales  by  viscous  heating.  This 
cascading  of  kinetic  energy  occurs  partially  within  a  wavenumber  space  (ranging  from 
only  a  few  centimeters  to  hundreds  of  meters)  referred  to  as  the  "inertial  sub-range".  It 
is  somewhere  within  this  inertial  sub-range  where  turbulent  scales  of  the  multiple  of 
one-half  the  transmitted  wavelength  are  encountered.  When  an  electromagnetic  wave 
front  encounters  these  areas  of  variability,  the  wave  front  is  refracted  slightly  and  a  small 
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amount  of  energy  is  scattered  in  all  directions,  as  represented  by  the  dotted  wave  fronts 
emanating  from  point  'AA'  within  turbulent  layer  'A'  in  Fig.  4.1  Very  high  above  the 
troposphere  the  density  of  air  is  such  that  it  cannot  support  small  scale  turbulence, 
hence  the  first  theoretical  limitation  on  the  profiler  is  a  maximum  altitude  that  can  sup¬ 
port  turbulence  on  the  scale  of  37  cm-a  height  which  varies  but  is  estimated  to  be  be¬ 
tween  18-20  km  for  the  West  coast. 

The  degree  of  turbulence  and  the  variations  in  refractivity  index  over  short  distances 
affect  the  energy  backscattered  (the  power  returned).  The  refractive  index  for  the 
troposphere  at  UHF  frequencies  (Balsley  and  Cage,  1980)  is 

/i=I-f77.6x  10*^  Y+ 3.73x10“'-^,  (1) 

where  e  is  the  partial  pressure  of  water  vapor  in  mb  (humidity),  p  is  the  atmospheric 
pressure  in  mb,  and  T  is  the  absolute  temperature  in  K.  A  refractive  index  structure 
function  Z),(r)  described  by  Gossard  and  Strauch  (1983)  is  used  to  quantify  these  vari¬ 
ations  over  short  distances  as 

(2) 

where  x  and  j?  -h  v  are  the  vectors  originating  at  the  profiler  antenna  and  terminating  at 
the  centers  of  two  air  parcels  (depicted  as  circles)  separated  by  a  distance  and  direction 
f ,  represented  in  Fig.  5.  Within  the  inertial  sub-range  of  homogeneous,  isotropic  tur¬ 
bulence,  Gossard  and  Strauch  (1983)  state  that 

Zl,(7)  oc  (3) 

To  make  an  equality,  a  standard  parameter  Q  ,  or  the  "refractive  index  structure  pa¬ 
rameter",  was  conceived.  Q  is  a  constant  of  proportionality  describing  the  refractive 
index  differences  between  two  parcels  of  air  situated  adjacent  to  one  another;  so  now 
Eq.  (3)  becomes 

D,(r)=Cy<^.  (4) 


1  There  is  evidence  suggested  by  Doviak  and  Zmic  (1984)  and  Waterman  et  al.  (1985)  that  at 
certain  wavelengths,  the  backscattered  energy  is  not  isotropic,  but  rather  anistropic.  There  is  on¬ 
going  research  to  resolve  this  phenomenon. 
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Fig.  4.  Illustration  of  the  Generation,  Beam  Pattern,  and  Scattering  of  Profiler 
Emiltcd  Signals 

Q  Is  dependent  on  the  temperature  and  humidity,  as  well  as  the  pressure,  as  given  by 
the  relationship  from  Peterson  (1988) 


C^=  79  X  10"^[-^][Cf  +  .1 13C7-2+  3.2  ^ 


Fig.  5.  Illustration  of  Refractive  Index  Structure  Sampled  by  the  Profiler  for  a 
Singic-Scattcrcr 

where  Q-  is  the  “temperature  structure  function”,  Cq  is  the  “humidity  structure  function”, 
and  Cj(,  is  a  combined  structure  function.  In  the  inertial  sub-range,  the  radar 
"rcllcctivity”  (tj)  is  proportional  to  Q,  and  is  given  by  Battan  (1973) 

1/  =  0.38CJr/l"‘'^  (6) 

where  A  is  the  radar  operating  wavelength. 

Radar  reHcctivity  can  now  be  expressed  as  returned  power  (P,),  derived  as  Uq.  (B-22) 
in  Appendix  B 


=  (7) 

R 

where  K  is  a  "radar  constant"  (i.e.,  radar  cfTicicncy  and  mathematical  constant),  P,  is  the 
power  transmitted,  t  is  the  pulse  duration,  A,  is  the  clTcctive  antenna  area,  R  is  the  dis¬ 
tance  from  the  antenna  to  the  turbulent  eddy,  and  tj  is  the  radar  reflectivity  from  Cq.  (6). 
R  is  determined  by  the  formula  R  =  ct/2,  i.e.,  the  lime  it  takes  for  backscattcrcd  energy 
from  the  respective  sampling  heights  to  return  to  the  antenna;  c  is  the  propagation  speed 
of  energy  at  radio  frequencies  and  t  is  the  time  of  transmission  subtracted  from  the  time 
of  receipt. 
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B.  PROFILER  DATA  PROCESSING 

The  wind  speed  and  direction  can  be  derived  trigonometrically  from  the  radial  ve¬ 
locity  of  the  profiler's  individual  transmitted  beams,  of  which  there  are  three-one  verti¬ 
cal  and  two  off-vertical.  Figure  4  graphically  depicts  how  these  beams  are  formed;  the 
off- vertical  beams  are  phased  by  delaying  the  transmitted  RF  pulse  to  certain  portions 
of  the  array.  'The  radial  velocity  is  derived  from  the  Doppler  frequency  shift  created  by 
turbulent  eddies  traveling  along  with  the  mean  wind.  In  the  presence  of  turbulence  or 
other  velocity  inhomogeneities  (refractive  index  fluctuations)  within  the  beam,  the  return 
signal  spectrum  may  vary  appreciably  from  one  pulse  to  the  next.  In  order  to  measure 
the  prevailing  flow  velocity  (by  averaging  out  microscale  fluctuations),  many  samples 
(about  1400)  are  composited  over  a  period  of  integration  (i.c.,  sample  period).  Even  if 
there  were  no  velocity  variations  within  the  beam,  the  returned  signal  would  be  at  an 
assortment  of  frequencies,  since  the  transmitted  signal  is  not  of  uniform  frequency.  (In 
addition  to  the  404  MHz  frequency,  the  transmitter  also  leaks  low-intensity  "tails"  at 
slightly  higher  and  lower  frequencies.  As  the  pulse  length  decreases  the  bandwidth  in¬ 
creases.) 

The  period  of  integration  is  58  seconds  in  the  low  mode  and  about  55  seconds  in  the 
high  mode.  The  result  at  the  end  of  the  integration  period  is  the  most  fundaniental 
profiler  quantity,  a  vertical  profile  of  36  wind  velocity  spectra.  A  representative  example 
of  the  velocity-power  spectra  is  given  in  Fig,  6.  There  are  six  modes  of  data:  two  off- 
vertical  and  one  vertical  mode  for  each  the  high  power  and  low  power  configurations. 
In  order  to  sample  the  entire  troposphere  and  lower  stratosphere,  low  and  high  mode 
profiles  are  measured  sequentially,  requiring  a  total  time  of  approximately  six  minutes. 
Thus,  10  sets  of  tropospheric  winds  can  oe  obtained  p3r  hour. 

To  get  data,  a  radio  wavefront  or  RF  pulse  must  first  be  transmitted.  After  the  en¬ 
ergy  is  radiated  by  the  antenna  system,  an  electronic  circulator  switch  is  thrown  to  re¬ 
connect  the  antenna  to  the  receiving  equipment.  There  are  100-150  microseconds 
between  each  pulse  and  it  is  during  this  interval  that  the  receiver  samples  the  returned 
signal  at  different  "gates"  or  vertical  range  bands  as  illustrated  at  the  top  of  Fig.  7. 
Then,  for  each  gate,  several  thousand  samples  are  gathered  and  time  domain  averaged 
by  the  processor,  as  shown  in  the  second  and  third  functional  steps  in  Fig.  7. 

Time  domain  averaging  (TDA)  consists  of  taking  the  mean  of  several  consecutive  sam¬ 
ples  in  time.  Figure  8  functionally  illustrates  how  the  received  energy  is  broken  into 
quadrature  and  in-phase  components,  and  the  subsequent  time-averaging  of  each  com¬ 
ponent.  This  time  domain  averaging  is  applied  to  the  power  spectra  of  the  36  available 
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Fig.  6.  Typical  Vertical  Profile  of  36  Range  Gate  Frequency-Power  Spectra:  (As 
displayed  on  the  Keyboard  Display  Unit  (KDU).)  In  the  upper-left  hand 
corner  of  each  gate  is  the  center  height  of  the  sampled  volume  in  meters. 
Beneath  the  height  parameter  is  the  resolved  radial  velocity  in  m  s'\  also 
indicated  by  the  vertical  bar  centered  within  the  dominant  peak.  In  the 
upper-right  of  each  spectrum  is  listed  the  signal-to-noise  and  below  this 
parameter  is  the  ma.ximum  pou.r  encountered  within  the  spectrum,  both 
in  dli. 


gates,  hence  a  TDA  matrix  can  be  formed  as  illustrated  in  Fig.  8,  which  consists  of  256 
samples  in  each  of  the  36  gatcs"256  samples  arc  needed  for  the  fast  fourlcr  transform 
(FFT)  calculations  later.  The  entire  time  the  processor  takes  to  collect  a  full  matrix  for 
all  range  gates  is  approximately  three  seconds  in  the  low  mode. 

'fo  this  point  in  the  data  processing,  a  typical  power  spectrum  in  relation  to  time  is 
available  for  each  gate.  But  what  is  desired  is  the  dominant  frequency,  therefore  fre¬ 
quency  domain  averaging  or  spectral  averaging  subsequently  ensues.  It  consists  of  cal¬ 
culating  a  mean  spectrum  from  several  consecutive  spectra  after  the  FFf  computations 
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Schematic  Diagram  of  the  Primary  Signal  Processing  Steps  for  the  Profiler 
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Fig.  8.  Funtional  Diagram  of  the  Time  Domain  Averaging  (TDA)  Process 


have  been  performed.  Now  the  time  required  to  complete  a  spectrum  in  each  of  the  six 
modes  can  be  determined: 


Time _per_spectrum  =  PRPx  TDA  x  256, 


where  PRP  is  the  pulse  repetition  period  of  the  radar  and  TDA  is  the  number  of  time 
domain  averaging  matrices.  Finally,  the  number  of  spectra  utilized  in  the  averaging 
process  of  each  mode  is  adjusted,  so  the  total  time  per  spectrum  is  approximately  60 
seconds,  and,  consequently,  the  total  lime  per  profile  is  then  six  minutes,  which  matches 
the  National  Oceanic  and  Atmospheric  Administration  reporting  requirements  for  their 
future  profiling  network.  The  number  of  spectra  utilized  for  a  single  profile  during  the 
first  month  of  operation  is  listed  in  Table  1. 

At  this  point,  a  typical  power  spectra  in  the  frequency  domain  has  been  processed 
for  each  mode  at  each  gate  from  which  wind  velocities  must  be  determined.  Briefiy,  a 
wind  profiler  is  a  Doppler  radar  that  measures  the  shift  in  frequency  created  by  moving 
targets.  In  the  case  of  the  wind  profiler,  the  targets  are  turbulent  eddies  moving  with  the 
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Table  1.  PROFILER  PARAMETERS 


Mode 

High 

Low 

No.  of 
TDA  Ma¬ 
trices 

40 

120 

PRP 

162ps 

100/rs 

Tinte  per 
Spectrum 

1.65s 

3.072s 

No.  of 
Spectra 

34 

19 

mean  wind.  Tims  by  determining  the  Doppler  sliift  of  the  spectral  moments,  the  wind 
toward  and  away  from  the  radar  can  be  determined.  To  accomplish  this,  a 
manufacturer's-supplicd  algorithm  computes  the  Doppler  spectral  moments  of  the 
spectral  peak  a.ssociated  with  a  local  maximum  in  the  frequency  spectrum.  Inputs  to  this 
algorithm  are  the  velocity  (or  frequency)  versus  power  spectrum  (with  noise  and  ground 
clutter  removed)  and  the  index  of  the  local  maximum  determined  to  be  signal.  Outputs 
of  this  algorithm  arc  the  three  spectral  moments: 

1)  Signal  Power  (estimated  total  power  of  the  received  signal), 

2)  Mean  radial  velocity  (first  moment  of  the  spectral  peak),  and 

3)  Radial  velocity  variance  (second  initial  moment,  or  width,  of  the  spectral 
peak). 

I  hi.s  routine  is  specified  for  each  mode  and  for  each  gate.  The  Doppler  shift  of  the  peak 
(f„)  can  be  converted  into  a  radial  velocity  component  of  the  wind,  by  rc-afranging  the 
following  basic  equation  from  Battan  (1979): 


'fhe  central  frequency  (as  shown  in  Fig.  7)  corresponds  to  the  transmitted  frequency. 
If  a  "target"  is  moving  toward  the  radar,  the  frequencies  will  be  incrcqsed;  conversely,  if 
it  is  moving  aw'ay,  the  frequency  will  decrease.  The  Doppler  return  and  the  width  of  the 
peak  at  each  range  gate  arc  measures  of  the  degree  of  turbulence  in  the  sample  volume. 
'I  he  second  peak  within  the  sample  spectra  shown  in  Fig.  7  is  an  indication  of  detected 
precipitation.  (Precipitation  would  be  verified  by  on-site  observations.) 
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Radial  velccity  components  are  next  decomposed  into  easterly  and  northerly  com¬ 
ponents  of  the  Cartesian  coordinate  system.  This  is  accomplished  by  the  vertical  and 
two  off-vertical  beams.  The  vertical  beam  pattern  measures  the  actual  w-component  of 
the  wind.  With  the  w-component  known,  the  u-  and  v-components  can  then  be  derived 
by  trigonometry,  as  shown  in  Appendix  C.  Within  Appendix  C,  and  6j  (the  angles 
of  the  offset  beams)  are  most  commonly  set  at  15*  from  vertical  and  tf)  (the  angle  in  the 
horizontal  plane  between  the  two  off-vertical  beams)  is  set  at  90®  for  case  of  trigono¬ 
metric  calculations. 

Finally,  the  three  moments  are  consensus-averaged,  based  upon  data  collected  for 
a  one-hour  averaging  period.  The  consensus  technique  is  detailed  in  Appendix  D.  As 
suggested  by  the  name,  not  all  the  data  are  averaged,  rather  only  those  that  mutually 
agree  within  a  threshold.  A  2  m  jt*  window,  or  threshold,  is  used  by  the  manufacturer 
for  the  radial  velocities.  If  the  consensus  group  does  not  have  at  least  four  members, 
then  the  average  value  is  reported  as  missing.  This  subset  is  then  used  to  calculate  the 
average  values  of  all  three  moments. 

C.  PROCESSING  OF  RAWINSONDE  DATA 

The  rawnsonde  signal  is  sampled  by  the  ground  system  every  ten  seconds.  The  wind 
data  from  24  consecutive  ten-second  samples  arc  then  averaged-this  provides  smoothing 
of  the  Initial  data  to  remove  noise  (see  Fig.  9).  Consequently,  all  sharp  features  are  ef¬ 
fectively  removed.  A  running  average  is  subsequently  maintained  wnth  data  output  every 
40  m,  based  on  a  four  mr'  balloon  ascent  rate.  A  detailed  presentation  of  the  scheme 
for  deriving  rawinsondc  data  carmit  be  given  within  this  study  because  the  material  is 
proprietary  with  Vaisala  Corporation.  (Winds  in  the  lowest  200  m  are  handled  in  a  quite 
different  fashion,  due  to  the  obvious  lack  of  data  required  for  a  running  average;  this 
different  technique  is  proprietary  information  but  is  not  important  for  this  study  since 
there  are  no  profiler  wind  values  at  these  levels.) 

The  rawinsonde  data  then  represent  approximately  a  1000  m  (actually  960m)  thick 
volume  (with  40  m  gate  spacing  between  subsequent  data),  whereas  the  wind  profiler 
has  250  m  resolution,  or  interval  spacing,  in  the  low  mode  (with  250  m  gate  spacing)  and 
1000  m  resolution  in  the  high  mode  (with  500  m  gate  spacing),  hence  different  values  as 
illustrated  in  Fig.  9.  Therefore,  a  typical  height  of  rawinsonde  measurements  centered 
within  the  sampling  profiler  gate  was  required  for  comparison.  U  and  V-velocity 
rawinsonde  components  were  subsequently  averaged  through  a  column  centered  around 
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the  desired  profiler  gate  height.  Likewise,  rawinsonde  temperature  and  humidity  were 
linearly  averaged  with  respect  to  height  around  the  profiler  gate  heights. 

D.  KNOWN  SOURCES  OF  PROFILER  ERROR 

There  are  many  potential  sources  of  noise  in  the  profiler  radial  velocities.  Of  most 
significance  for  this  comparison  study,  the  profiler  and  rawinsonde  sample  at  different 
locations.  The  wind  profiler  provides  simultaneous  measurements  over  all  heights  at  a 
fixed  location.  It  samples  large  volumes  of  atmosphere  (up  to  1  km  deep)  and  assigns 
an  integrated  value  to  the  center  height  of  the  sampled  column.  The  rawinsonde,  on  the 
other  hand,  records  measurements  at  different,  discreet  heights  and  at  different  times  (a 
typical  ascent  will  take  up  to  1-1/2  hours  and  may  cover  many  kilometers  lateraUy).  As 
a  result,  the  rawinsonde  and  the  wind  profiler  often  obtain  measurements  at  widely  sep¬ 
arated  locations  where  the  winds  may  be  different.  Also,  the  winds  over  the  fixed  profiler 
location  are  not  uniform  during  the  rawinsonde's  1-1/2  hour  ascent;  therefore,  variability 
of  the  atmosphere  can  make  the  interpretation  of  any  comparisons  difficult  at  times. 
Furthermore,  since  these  instruments  sample  in  a  functionally  different  manner,  com¬ 
parisons  of  their  measurements  cannot  provide  unambiguous  information  on  the  accu¬ 
racy  of  either  instrument. 

Another  large  source  of  error  is  electromagnetic  noise  present  within  the  pulse  vol¬ 
ume  as  a  result  of  nearby  power  lines  or  machinery,  as  well  as  other  contributions  fror.i 
side  lobes  (power  “leaking"  from  the  antenna  at  angles  other  than  the  desired  angle  >. 
The  Federal  Aviation  Administration  began  operating  a  traffic  control  radar  within  60 
m  of  the  Ft  Ord  wind  profiler  site  in  April  of  1990.  A'ihough  the  operating  frequency 
docs  not  appear  to  affect  profiler  operations,  the  movement  of  the  antenna  provides  a 
strong  enough  Doppler  shift  of  the  profiler  transmitted  frequency  to  provide  a  significant 
return  in  the  side  lobes  of  the  receiver  entenna  and  contaminate  the  lowest  gate.  There 
appear  to  be  many  forms  of  this  contamination  such  as  the  inclusion  of  a  sinusoidal 
wave  pattern  in  the  lowest  four  gates  as  represented  in  Fig.  10.  A  more  common  con¬ 
tamination  pattern  is  the  occurrence  of  narrow  peaks  at  approximately  1ms-'  (not 
shown). 

Noisy  wind  vectors  can  result  from  an  imbalance  of  the  individual  beam  sensitivities, 
since  one  beam  may  .detect  the  weak  atmospheric  turbulent  signal  while  the  other  may 
intermittently  resolve  the  background  noise  as  the  dominant  signal  source.  In  the  event 
that  the  atmosphere  is  quite  dry  and  the  associated  profiler  backscatter  w'eak,  there  is  a 
fairly  strong  persistent  background  noise  frequently  present  and  the  “insensitive"  beam 
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Fig.  9.  Functional  Ran  insonde  and  Profiler  Gate  Spacing 


may  'lock  in  on"  the  background  frequency.  The  result  is  a  high  consensus  nuntber  but 
an  unrealistic  wind. 

Closely  related  to  this  latter  error  is  the  susceptibility  of  each  beam  to  be  less  sensi¬ 
tive  to  cross-winds,  or  those  prevailing  winds  which  cross  perpendicular  to  the  direction 
which  the  beam  is  formed.  Hence,  the  accuracy  of  determining  these  cross-wind  veloci¬ 
ties  is  reduced.  There  is  consistent  evidence  of  this  bias  throughout  the  course  of  this 
study,  and  analysis  of  a  larger  data  set  would  probably  define  the  bias  discreetly. 

When  the  radar  switches  from  transmit  to  receive  mode,  there  is  electronic  noise  in 
the  system  that  takes  a  few  microseconds  to  dissipate.  Any  signals  returned  by  the  at¬ 
mosphere  during  this  brief  period  are  masked  by  noise.  Because  the  first  signals  to  re¬ 
turn  come  from  low  in  the  atmosphere,  the  lowest  sampling  gate  (446  m)  is  most 
susceptible  to  this  problem.  No  evidence  of  this  elTcct  was  discovered  through  the 
course  of  this  study. 

The  profiler  operating  frequency  and  sampling  interval  determine  the  Nyquist  ve¬ 
locity,  the  maximum  radial  velocity  that  can  be  measured.  The  well-known  Nyquist 
sampling  theorem  states  that  the  highest  frequency  to  be  resolved  in  the  signal  must  be 
sampled  at  least  two  times  per  period.  If  the  signal  contains  energy  at  frequencies  higher 
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Fig.  10.  Sample  of  FA  A  Traffic  Control  Radar  Contamination  of  NFS 
Profiler:  Evident  in  the  four  lowest  gates.  (See  Fig.  6  for  format  de¬ 
scription.) 


than  the  Nyquist  frequency,  the  Doppler  signal  can  be  aliased  by  higher  frequencies, 
'fhis  means  that  the  energy  is  "folded  back"  to  frequencies  less  than  the  Nyquist  fre¬ 
quency  (McGlllcn  and  Cooper,  1974).  For  the  NFS  profiler  the  Nyquist  velocity  is 
about  15  ms~'  (30  kt)  in  the  low  mode  and  60  ms*‘  (120  kts)  in  the  high  mode.  Radial 
velocities  greater  than  the  Nyquist  velocity  arc  "folded"  as  illustrated  in  Fig.  1 1  at  1 1697 
m.  Although  velocity  folding  occurred  in  two  spectra  within  Fig.  1 1,  fortunately  it  was 
not  folded  into  the  spectral  region  "sensitive"  to  the  algorithm.  There  were  no  known 
occurrences  of  this  type  error  during  the  study  period. 

Contributions  from  persistent  vertical  velocity,  such  as  in  the  vicinity  of  sloping 
terrain  or  the  closed  cells  of  the  land/sea  breeze  circulation,  may  cause  appreciable  con¬ 
tamination  of  the  "horizontal"  velocity  measured  by  the  off-zenith  beams,  especially 
when  the  horizontal  wind  speed  is  low.  Since  the  "horizontal"  beams  are  15®  from 
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Fig.  II.  Sample  of  Velocity  Folding  in  the  NFS  ProFilen  Seen  at  the  11697  m 
level.  (See  Fig.  6  for  format  description.) 

zenith,  75  percent  of  the  oflT-zenith  radial  velocity  is  sensed  as  vertical  velocity.  Thus  the 
horizontal  velocity  required  to  produce  a  I  mj~'  radial  component  is  about  4  nts’K 
Synoptic-scale  vertical  velocities  (~  0.1  ms*')  appear  not  to  pose  a  problem,  but 
smaller-scale  systems  which  may  contain  vertical  velocities  as  great  as  several  meters  per 
second  arc  clearly  a  potential  source  of  horizontal  beam  contamination. 

Turbulence  accompanied  by  significant  backscatter  above  a  height  which  the  profiler 
can  "see"  introduces  aliasing  into  the  data.  The  ma.ximum  unambiguous  radar  height  is 
given  by  Battan  (1974)  as 


R, 


ciPRP) 


mu 


where  c  is  the  propagation  speed  of  radio  energy.  PRP  as  set  in  the  low  mode  of  the 
NPS  profiler  is  100  fis,  which  sets  in  the  low  mode  at  15,000  m.  (/?„„  in  the  high 
mode  is  24,300  m.)  Therefore,  if  the  transmitted  signal  is  returned  from  a  source  higher 
than  15,000  m,  the  radar  processor  while  in  the  low  mode  will  sec  this  signal  as  a  possible 
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early  return  in  the  lower  gates  of  the  next  sampled  pulse.  Consequently,  ionospheric  and 
cosmic  noise  can  also  b,e  problems.  In  this  case  the  radar  will  receive  noise  in  any  of  the 
gates,  when  realistically  the  signal  evolves  quite  a  bit  higher  beyond  that  which  the  radar 
can  measure. 

Precipitation  and  other  hydrometeors  can  also  scatter  energy  back  towards  the  an¬ 
tenna.  Typical  meteorological  radars  operate  at  wavelengths  on  the  order  of  5  to  10  cm, 
whereas  wind  profilers  operate  typically  near  wavelengths  of  74  cm.  Hydrometeors  are 
some  10*  times  less  reflective  to  wind  profilers  compared  to  conventional  weather  radars. 
This  is  indicated  by  the  reflectivity  equation  given  by  Gossard  and  Strauch  (1983)  for 
single  particle  scattering: 


{const)t^ 


where  d  is  the  particle  diameter.  Precipitation,  when  detected,  can  be  a  problem  in  some 
circumstances.  Clearly,  if  there  is  interest  in  the  vertical  motion  of  the  air,  precipitation 
presents  a  problem,  as  the  vertical  velocity  measured  is  that  of  the  precipitation,  and  not 
the  air.  In  some  cases,  spectra  from  the  vertical  beam  will  show  a  double  peak:  one 
corresponding  to  the  clear  air  and  the  other  precipitation.  This  is  illustrated  in  Fig.  12. 

Rawinsondc  operation  colocated  with  the  profiler  creates  a  false  wind  velocity  value. 
Figure  13  presents  two  spectra  from  the  Keyboard  Display  Unit.  Balloon  interference 
is  highlighted  at  approximately  6000  m  in  the  top  frame  and  11000  m  in  the  bottom 
frame.  The  time  diflerence  between  these  two  frames  is  25  minutes.  Again,  assuming  a 
A  rnsr'  balloon  ascent  rate,  this  is  precisely  the  height  the  rawinsonde  should  have  at¬ 
tained.  The  detection  of  the  balloon  appears  to  be  due  to  the  different  index  of  refrac¬ 
tion  value  characteristic  of  the  gas  within  the  balloon  vice  the  actual  reflection  from  the 
non-mctallic  surface  of  the  balloon.  The  radio  signal  of  402  MHz  utilized  by  the 
rawinsonde  for  passing  data  to  the  ground  system  does  not  appear  to  interfere  with 
profiler  operation. 

It  must  also  be  recognized  that  the  vertical  beam  and  the  two  off-zenith  beams  do 
not  sample  the  same  volume  of  air,  in  either  space  or  time.  There  may  be  up  to  a  3.2 
km  difference  at  an  altitude  of  12  km.  This  is  especially  troublesome  in  the  presence  of 
small-scale  mesoscale  disturbances. 
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Fig.  12.  Example  of  a  Double  Peak  in  NPS  Profiler  Spectra  causediiy  Rain;  ev¬ 
ident  clearly  at  the  914  and  1214  m  levels,  as  well  as  other  disparities  ex¬ 
tending  up  to  4000  m.  (Sec  Fig,  6  for  format  description.) 
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in.  DATA  ANALYSIS 


A.  COMPARISON  OF  MANUFACTURER'S  CONSENSUS  TO  RAWINSONDE 
DATA 

Prior  to  analyzing  dynamic,  kinematic  and  thermodynamic  relationships  derived 
from  the  NPS  profiler  winds,  it  is  necessaiy  to  first  determine  the  accuracy  of  the  profiler 
wind  data.  This  information  provides  the  NPS  Marine  Atmospheric  Measurements  Lab 
a  benchmark  of  how  the  NPS  profiler  performs  in  relation  to  other  profilers  as  well  as 
verifying  requirements  within  the  profiler  procurement  contract.  Both  the  wind  profiler 
and  rawinsonde  system  provide  vertical  profiles  of  wind  measurements;  therefore,  like 
many  previous  studies,  the  accuracy  of  the  profiler  will  be  gaged  by  the  rawinsonde  data 
taken  as  truth  data.  A  statistical  comparison  between  profiler  and  rawinsonde  wind  data 
was  made  and  used  two  approaches:  (I)  each  horizontal  wind  component  was  correlated 
between  the  profiler  and  rawinsonde  as  determined  by  a  correlation  coefficient,  and  (2) 
mean  and  standard  deviations  for  the  difference  in  each  horizontal  wind  component  as 
measured  by  the  two  instruments  was  calculated.  Because  such  a  small  rawinsonde  data 
set  was  available  during  the  first  month  of  operation,  a  more  detailed  look  into  why 
potentially  "bad"  profiler  data  was  not  in  agreement  with  the  rawinsonde  data  was  done 
manually,  providing,  a  third,  more  subjective  approach. 

Nine  individual  profiles  are  compared  to  assess  the  quality  of  the  manufacturer's 
consensus  winds  measured  by  the  profiler.  Figure  14  is  a  representative  comparison  of 
the  vertical  profiles  of  wind  speed  and  direction  measurements  for  20  April  1990  made 
by  the  two  instruments.  The  associated  complete  rawinsonde  sounding  is  presented  in 
Fig.  15.  Additional  comparisons  for  the  remaining  eight  soundings  (Figs.  E-I  through 
E-16)  as  w'ell  as  the  specific  positions  of  the  five  at-sea  rawinsonde  launchings  (Fig. 
E-17)  arc  presented  in  Appendix  E.  In  each  case  the  rawinsonde  profile  is  flanked  by  its 
associated  Tycho  hourly-consensus  profiles,  approximately  2  hours  on  either  side  of  the 
rawinsonde  launch  lime.  Subjectively,  the  profiles  correlate  very  well,  at  least  where 
consensus  in  the  radar  data  is  achieved.  The  wind  direction  correlates  very  well.  Below 
5  km  the  wind  speeds  are  generally  within  5  kts  of  each  other,  while  above  5  km  the  wind 
speeds  are  accurate  to  within  10-15  kts.  The  departure  between  individual  profiler  and 
rawinsonde  wind  barbs  in  the  lowest  2000  m  seen  from  the  remaining  launch  times  in 
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Appendix  E  could  potentially  be  due  to  the  discrepancy  in  instrument  locations  as  well 
as  the  turbulence  detected  by  the  radar  profiler  in  the  presence  of  the  marine  inversion. 

A  comparison  (rawinsonde  vs  profiler)  of  the  individual  horizontal  wind  components 
versus  height  are  illustrated  in  Fig.  16  and  17  for  a  representative  low  mode  profile  from 
0000  UTC  20  April  1990  and  Fig.  18  and  19  for  a  typical  high  mode  profile  from  0700 
UTC  03  May  1990.  There  is  generally  very  good  agreement  below  5-6  km  in  the  low 
mode  and  below  8-9  km  in  the  high  mode.  There  is  a  notable  bias  in  the  v-velocity 
component  in  the  high  mode,  which  appears  to  be  due  to  the  orientation  of  the  prevail¬ 
ing  winds  in  relation  to  the  profiler  antenna,  which  in  this  case  were  from  the  north. 

In  the  low  mode  the  maximum  height  obtained  appears  to  correlate  with  the  fall-off 
in  humidity  (evident  in  Fig.  15  where  the  atmosphere  becomes  dry  above  4000  m);  but 
is  never  greater  than  approximately  8  km  even  when  humid  layers  are  present  aloft.  This 
is  due  to  the  weakened  signal  from  geometric  spreading.  (A  humid  layer  aloft  is  defined 
in  this  paper  as  a  layer  with  greater  than  40%  relative  humidity  above  4000  m.)  This 
limitation  is  expected  since  a  .74  cm  radio  signal  is  four  times  more  sensitive  to  humidity 
fluctuations  associated  with  turbulence  than  to  temperature  fluctuations,  at  least  in  the 
lower  atmosphere  (Cage  1990).  This  is  illustrated  in  Fig.  20,  by  a  deep  humid  layer 
which  exists  to  5  km.  This  is  reflected  in  both  the  U  and  V-velocity  component  winds 
in  Fig.  18  and  Fig.  19. 

As  additional  proof  of  the  dominance  of  humidity  in  the  lower  atmosphere,  a  nearby 
sounding  was  compared  with  the  raw  spectra.  Fig,  20  is  the  1200  UTC  06  April  1990 
Oakland  sounding,  10  hours  prior  to  the  low-power  mode  profiler  spectra  in  Fig.  21 
(0000  UTC  07  April  1990  was  not  available).  Typical  conditions  for  this  period  of  study, 
as  seen  in  Fig.  21,  are  first  a  layer  of  very  strong  returns  between  500-1400  m,  the  ma¬ 
rine  thermal  inversion;  which  is  followed  by  a  rapid  fall-off  of  signal  return,  until  ap¬ 
proximately  3000  m  is  reached  where  the  signal  levels  begin  to  increase  due  to  humidity. 
This  corresponds  very  well  with  the  moisture  layers  in  the  sounding.  This  is  then  fol¬ 
lowed  by  a  fall-off  between  4000-5000  m  and  a  local  maximum  of  -2  dB  S/N  reached  at 
5000-5500  m.  In  summary,  a  very  good  subjective  correlation  between  the  rawinsonde 
sounding  and  profiler  spectra  is  evident. 

The  maximum  height  attained  in  the  high  mode  is  about  12  km  (significantly  less 
than  the  theoretical  24  km  which  the  electronics  allow),  again  due  to  weakened  signal 
return  from  geometrical  spreading.  But,  what  is  additionally  seen  in  the  U  and  V  com¬ 
ponents  is  a  2-3  km  gap  in  consensus  data  just  below  the  tropopause  (see  Fig.  18  and 
Fig.  19).  The  tropopause  is  identified  by  a  constant  temperature  profile  seen  in  the 
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rawinsonde  data.  It  appears  at  this  time,  that  this  gap  is  centered  around  the  maximum 
of  the  "jet  stream",  or  jet  stream  core  if  overhead.  This  is  the  most  consistent  feature 
seen  in  the  data  above  8  km.  Apparently,  the  wind  shear  throughout  this  region  is 
minimal,  thus,  reduced  turbulence  in  a  relatively  dry  region  provides  little  signal  return. 
However,  above  the  jet  maximum,  winds  begin  to  decrease  and  radar  returns  are  de¬ 
tectable  above  noise.  This  too  is  seen  in  the  raw  spectra  in  Fig.  22,  where  signal 
strength  falli  off  and  then  increases  above  -10  dB  at  the  12197  m  level. 

To  quantify  the  data  accuracy,  a  statistical  comparison  was  made  by  computing  the 
correlation  cocITicient,  calculating  the  standard  deviation,  and  graphically  portraying  the 
data  "pairs"  on  scatter  plots.  The  comparisons  were  based  on  unedited  consensus  data 
from  the  manufacturer's  algorithm.  Figures  23  and  24  are  scatter  plots  of  component 
rawinsonde  measurements  versus  wind  profiler  measurements  for  the  nine  combined 
launchings.  Each  marker  corresponds  to  a  measurement  made  at  about  the  same  hour 
and  height  for  both  instruments.  If  the  two  measurements  are  equal,  then  the  marker 
will  lie  on  the  diagonal.  If  the  marker  is  positioned  off  this  diagonal  the  two  measure¬ 
ments  arc  different.  Consistent  markers  to  one  side  or  the  other  indicate  a  bias. 

In  the  low  mode,  the  density  of  data  around  the  scatter  plot  diagonal  appears  good 
below  7-8  (15-17kts).  Winds  above  this  value  are  generally  associated  with  higher 

heights;  and  as  Fig.  23b  illustrates,  the  "bad"  values  (the  values  that  do  not  fall  closely 
along  the  diagonal)  generally  lie  above  6  km  as  represented  by  the  markers  O  and  A. 
Similarly,  in  the  high  mode  data  depicted  in  Fig.  24b  there  appears  to  be  an  appreciable 
density  of  data  along  the  diagonal,  uith  any  "bad"  data  related  to  speeds  of  8  m  j-* 
( 1 7kts)  above  8  km.  The  depiction  of  "good"  data  at  higher  speeds  indicate  that  height 
is  the  primary  contributor  to  the  presence  of  "bad"  data.  There  is  a  bias  tendency  in 
v-velocity  components  which  is  more  noticeable  in  the  low  mode  data.  This  cannot  be 
explained  but  may  be  due  to  such  a  small  data  set,  or  due  to  the  orientation  of  the  pre¬ 
vailing  winds  relative  to  the  profiler  fixed  beams. 

Table  2  presents  the  statistical  results  broken  out  by  atmospheric  layers.  The  values 
listed  under  the  column  "Mean",  represent  the  mean  error  away  from  the  mean  wind 
component,  while  "Standard  Deviation",  is  that  value  taken  around  the  mean  wind 
component,  within  one  standard  deviation  or  67%  of  the  data.  The  "Correlation  Coef¬ 
ficient"  values  were  calculated  by  the  Pearson  correlation  coefficient  relating  profiler  to 
rawinsonde  data.  'The  guidelines  this  study  has  set  for  "good"  data  is  a  correlation  co¬ 
efficient  above  .93  and  a  standard  deviation  less  than  2.5  m  s~'  ,  which  are  taken  from 
the  findings  by  Weber  and  Wuertz  (1990).  However,  the  performance  of  the  radar  in 
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relation  to  the  rawinsonde  system  never  quite  met  these  criteria.  The  table  reveals  good 
correlation  in  the  low  mode  below  6  km  in  both  components.  Good  correlation  in  the 
high  mode  is  seen  below  8  km.  There  is  poor  correlation  in  the  lowest  2-3  km  under  the 
presence  of  the  inversion  layer,  i.e.,  there  is  turbulence  created  by  the  wind  shear,  the 
inversion  jet,  and/or  convective  thermal  cells.  The  statistical  results  appear  to  reinforce 
the  conclusions  made  earlier  by  the  subjective  analysis. 

An  accuracy  determination  would  not  be  complete  without  the  inclusion  of  the  site 
survey,  which  establishes  exactly  where  the  antenna  is  oriented.  For  consideration,  the 
elevation  accuracy  of  the  antenna  base  height  is  ±  1  m,  while  the  azimuthal  accuracy  is 
±  3'  of  arc.  The  results  and  the  survey  site  layout  are  presented  in  Appendix  F.  The 
accuracy  to  which  the  profiler  array  was  established  is  considered  negligible  for  this 
study,  but  is  presented  for  completeness  and  for  the  record  to  accommodate  future 
studies. 

The  NFS  wind  profiler  appears  from  the  statistical  analysis  to  be  accurate  to  within 
only  3-4  m  r-'  in  the  lower  atmosphere  (or  more  specifically  within  the  moist  layer  below 
6  km).  This  error  is  slightly  greater  than  the  2.5  m  r-'  standard  deviation  computed  by 
Weber  and  Wuertz  (1990).  The  profiler  error  is  significantly  greater  than  the  1  ms"' 
accuracy  specification  listed  in  the  procurement  contract.  Meanwhile,  Strauch  et  al 
(1987)  found  a  standard  deviation  of  about  1.3  mj-'  in  clear  air  and  Wuertz  ei  al.  (1988) 
found  a  standard  deviation  of  about  2-4  m  r'  in  precipitation.  Thomson  and  Williams 
(1990)  reported  amazing  accuracies  in  profiler  data  of  ±l  msrK  The  studies  by  Weber 
and  Wuertz,  as  well  as  the  other  studies  presented,  analyzed  several  thousand  pairs  of 
rawinsonde  and  profiler  data,  while  this  study  only  was  able  to  look  at  a  few  hundred. 
A  larger  data  set  would  be  expected  to  reveal  a  lower  standard  deviation  based  on  the 
small  mean  error  values  listed  in  Table  2,  for  heights  below  6  km. 
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Sonde  High  Mode  Comparison  Data:  90  04  19  at  23  23  GMT 


o  K 
•<r  o 


O  o 

N  .. 


s 

C 

I 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

00 

(O 

o 

00 

u> 

N 

(sj3}9Uj)  apnijjiv 


Fig.  14.  Pronier*Ra>unsonde  Comparison  Winds  for  0100  UTC  20  April  1990 
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15.  Rawinsonde  Sounding  for  0100  UTC  20  April  1990 
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Fig.  16.  Typical  Comparison  Plot  of  Ran  insonde  vs  Profiler  U- Velocity  Component 
(Profiler  in  Low-Mode):  Rawinsondc  winds  are  indicated  by  the  o 
marker,  while  the  profiler-derived  winds  are  indicated  by  the  Hh  marker. 
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Fig.  17.  Typical  Comparison  Plot  of  Ran  insonde  vs  Profiler  V-Velocity  Component 
(Profiler  in  Lon-Mode):  Rawinsondc  winds  are  indicated  by  the  o 
marker,  while  the  profiler-derived  winds  are  indicated  by  the  +  marker. 
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Fig.  18.  Typical  Comparison  Plot  of  Rawiiisonde  vs  Profiler  U-Velocity  Component 
(Profiler  in  High-Mode):  Rawinsonde  winds  are  indicated  by  the  o 
marker,  while  the  profiler-derived  winds  are  indicated  by  the  +  marker. 
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Fig.  19.  Typical  Comparison  Plot  of  Rawinsonde  vs  Profiler  V-Velocity  Component 
(Profiler  In  High-Mode):  Rawinsonde  winds  arc  indicated  by  the  o 
marker,  while  the  profiler-derived  winds  are  indicated  by  the  +  marker. 
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Fig.  20.  OakJand  (72493)  Sounding  for  1200  UTC  06  April  1990 
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Fig.  21.  Example  of  Lo>v  Mode  Spectra:  There  is  increased  signal  strength  in  the 
5196  and  5446  m  levels  due  to  humiditj'. 
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Fig.  22.  Sample  Lotv  Mode  Spectra  Depicting  Jet  Stream  Structure:  Evident  is 
an  increase  of  signal  strength  at  the  12197  m  level. 
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Sonde  Speed  fnyi)  K*] 


ig.  23.  Scatter  Plot  of  Lo>v-Mode  Profiler  v&  Rawinsonde  Velocity:  Data  taken 
from  all  nine  rawinsonde  launches.  The  marker  x  represents  those  winds 
derived  below  2  km,  o  from  2  to  4  km,  +  from  4  to  6  km,  A  from  6  to  8 
km,  and  O  any  winds  above  8  km. 


Fig,  24.  Scatter  Plot  of  High-Mode  Profiler  vs  RaMinsonde  Velocity:  Data  taken 
from  all  nine  rawinsonde  launches.  The  marker  x  represents  those  winds 
derived  below  4  km,  o  from  4  to  6  km,  +  from  6  to  8  km,  A  from  8  to  10 
km,  and  O  any  winds  above  10  km. 
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Table  2.  STATISTICAL  ANALYSIS  OF  DIFFERENCES  BETWEEN 
RAWINSONDE  AND  PROFILER-DERIVED  VELOCITIES  FOR  ALL 
NINE  RAWINSONDE  LAUNCHES 


wm 

Layer 

mmSmM 

Standard 

Deviation 

Correlation 

Coefficient 

All  Heights 
Combined 

2.49 

10.31 

.490 

<2  km 

0.083 

3.4 

.791 

Low 

U  Velocity 

2-4  km 

-.828 

4.66 

.815 

4-6  km 

-.899 

8.4 

.783 

<6  km 

-.553 

5.78 

.793 

6-8  km 

-4.59 

14.09 

.285 

>  8  km 

-13.0 

18.68 

.077 

All  Heights 
Combined 

7.54 

14.31 

.415 

<2  km 

2.06 

4.04 

.634 

Low 

V  Velocity 

2-4  km 

1.67 

4.31 

.757 

4-6  km 

3.10 

8.01 

.714 

<6  km 

2.22 

5.65 

.740 

6-8  km 

18.98 

19.25 

0.01 

>8  km 

26.0 

20.34 

-.048 

I 

All  Heights 
Combined 

-4.94 

12.05 

.439 

<4  km 

-.996 

3.275  . 

.912 

High 

U  Velocity 

4-6  km 

-1.04 

3.091 

.970 

6-8  km 

0.627 

13.78 

.516 

‘ 

<8  km 

-.717 

6.44 

.805 

8-10  km 

-15.36 

21.43 

.063 

>  10  km 

-16.35 

14.69 

.647 

All  Heights 
Combined 

4.88 

12.83 

.361 

<4  km 

1.985 

3.817 

.742 

High 

V  Velocity 

4-6  km 

0.547 

3.916 

.914 

6-8  km 

4.29 

10.77 

.611 

<8  km 

1.99 

5.79 

.783 

8-10  km 

16.88 

29.61 

-.609 

>  10  km 

11.59 

18.96 

.217 
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B.  INVESTIGATION  OF  VARIABLES  AFFECTING  DATA  QUALITY 

Given  that  there  are  numerous  sources  of  error  in  the  profiler  operatioii  that  influ¬ 
ence  the  six-minute  profiles  and  hourly-consensus,  it  may  be  possible  to  eliminate  bad 
six-minute  data  in  order  to  subsequently  improve  overall  performance.  To  assess  the 
six-minute  data,  several  parameters  related  to  the  basic  moments  calculated  for  each  gate 
are  examined:  Signal-to-Noise  power  Ratio  (S/N),  beam  power  differential,  and  velocity 
variance  (or  width).  Although  one  of  the  most  effective  data  quality  control  techniques 
is  a  running  comparison,  this  study  will  attempt  to  identify  only  those  parameter 
thresholds  indigeneous  to  their  respective  sbt-minute  profiles,  in  the  hope  of  retaining 
as  many  mesoscale  atmospheric  features  as  possible.  This  study  also  proceeds  with  the 
presumption  that  by  improving  six-minute  profiles,  the  hourly  consensus  will  be  im¬ 
proved. 

In  order  to  place  any  credence  on  these  thresholds  it  was  necessary  to  analyze  the 
largest  data  base  possible.  Unfortunately  due  to  the  combination  of  power  outages, 
international  project  requirements,  and  equipment  maintenance  downtimes,  the  data 
base  is  limited.  Also,  in  order  to  compare  these  thresholds  the  direction  and  .speed  of  the 
prevailing  winds  would  have  to  be  the  same.  By  using  the  profiler  consensus  data,  a 
10-hour  period  on  02  May  1990  (0600-1600  UTC)  was  chosen  which  displayed  the  most 
uniformity  over  time  in  the  5000-7000  m  levels.  Based  on  findings  from  the  previous 
section,  this  is  the  highest  elevation  which  the  radar  can  feasibly  measure  reliably,  i.e., 
within  3-4  m  s-\  The  upper  atmosphere  (above  7  km)  was  assumed  to  be  less  variable 
than  the  lower  levels  (no  major  synoptic-scale  traveling  features  could  be  identified  in  the 
NMC  upper  air  analyses),  and  the  lower  atmospheric  layers  below  5  km  were  probably 
just  as  variable  as  at  any  other  time  period  which  could  be  chosen  (again,  there  were  no 
identifiable  wind  shifts  from  frontal  passages,  for  example,  which  could  be  identified). 
The  manufacturers  consensus  data  for  this  analysis  period  is  presented  in  Fig.  25,  and 
does  in  fact  indicate  very  uniform  conditions. 

The  sensitivity  of  the  receiver  electronics  is  expected  to  be  a  limiting  factor.  This  is 
expressed  as  a  ratio  of  returned  signal  strength  to  noise.  Signal-to-Noise  (S/N)  was  first 
calculated  for  each  of  the  profile  gates  from 

SfN  =  Power  —  [Noise  -b  10  log„(a))], 

where  Power  is  the  integrated  power  of  the  peak  signal,  Noise  is  the  signal  strength  of 
noise  under  the  peak,  and  to  is  the  bandwidth  to  compensate  for  the  noise  not  being  in¬ 
tegrated.  In  the  analysis,  (o  =  256,  for  256  different  sampling  points  in  the  frequency 
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Fig.  25.  Profiler  Hourly-consensus  Wind  Data  from  0605  to  1703  UTC  02  May 
1990:  This  consistent  period  represents  the  data  set  used  to  calculate 
specific  thresholds  from  the  spectral  moments. 
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spectrum.  Next,  the  U  and  V  wind  components  from  the  six-minute  files  were  decom¬ 
posed  into  the  original  radial  velocities.  Each  of  these  radial  velocities  are  plotted  versus 
their  associated  S/N  value  in  Fig.  26,  all  heights  included.  The  result  is  a  large  concen¬ 
tration  of  points  above  -19  dB,  and,  conversely,  a  very  broad  scatter  of  points  below  -19 
dB.  Statistically  the  standard  deviation  decreases  dramatically,  as  noted  in  Table  3.  The 
S/N  thresholds  in  Fig.  26  are  slightly  different  in  each  of  the  beams,  e.g.,  -18  dB  might 
be  used  as  a  value  instead  of  -19  dB-which  is  presumed  to  be  related  to  the  direction 
of  the  prevailing  wind,  out  of  the  north,  which  would  favor  the  X-beam.  Further  studies 
on  larger  data  sets  of  profiler  wind  data  could  better  define  this  threshold. 


Table  3.  S/N  THRESHOLD  STATISTICS 


Velocity 

Compo¬ 

nent 

<-21  dB 

-21to-20  dB 

-20to-19  dB 

-19to-18  dB 

>-18  dB 

X 

29.94 

23.86 

10.36 

8.84 

6.47 

Y 

41.72 

31.09 

24.45 

10.16 

10.20 

As  discussed  earlier,  the  profiler  resolves  horizontal  winds  by  sampling  large  volumes 
of  air  by  offset  beams  displaced  by  several  meters  from  one  another.  Therefore,  a 
threshold  is  needed  to  throw  out  data  where  supposedly  one  beam  will  lock  on  a  heli¬ 
copter,  airplane,  or  atmospheric  noise  and  possibly  even  ionospheric  noise.  Thus,  to 
relate  the  two  beams  a  power  difference  is  calculated,  simply  one  returned  power  value 
subtracted  from  the  other.  Figures  27  and  28  represent  the  power  difference  between  the 
X  and  Y-bcams,  for  above  5000  m  and  less  than  3000  m  respectively.  Above  5000  m  the 
beam  power  differential  docs  not  affect  the  mean  velocity  nor  the  standard  deviation  of 
the  data  points,  that  is  to  say,  the  density  of  markers  seem  to  be  equally  dispersed 
around  the  mean  wind  velocity  component.  This  is  presumably  because  the  signals  are 
relatively  so  much  weaker.  Below  3000  m,  however,  a  shift  in  the  mean  velocity  is  evi¬ 
dent  when  a  5  dB  differential  is  exceeded;  the  mean  velocity  and  standard  deviation  are 
plotted  out  separately  in  Fig.  29.  In  the  X-beam,  the  shift  in  mean  velocity  is  quite  ev¬ 
ident  above  5  dB,  but  not  so  evident  in  the  Y-bcam  (a  shift  of  only  1  m  at  6  dB).  The 
standard  deviation  throughout  the  plots,  represented  by  the  length  or  range  of  the 
shaded  bars,  show  no  apparent  trend.  Again,  the  X-beam  statistics  are  presumed  more 
prominent  due  to  the  orientation  to  the  prevailing  winds.  Subjectively,  by  eyeballing 
individual  wind  barb  plots  which  highlight  suspected  "bad"  winds,  the  5  dB  value  can 
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be  refined  to  6  dB,  i.e.,  more  good  data  in  relation  to  the  consensus  is  retained.  After 
all,  an  objective  of  this  study  is  to  make  the  six-minute  profiles  as  representative  and  as 
useful  as  possible. 

Finally,  the  velocity  variance  was  examined.  This  proved  the  most  difficult  thresh¬ 
old  to  determine.  The  premise  was  that  the  manufacturer's  algorithm  could  not  distin¬ 
guish  between  two  different  peak  velocities,  and,  consequently,  would  average  out  the 
two  peaks-efiectively  masking  the  true  wind  velocity  by  the  presence  of  another  signal, 
such  as  rain  or  helicopters.  This  is  illustrated  by  Fig.  30  at  914  to  1214  m  in  the  spec¬ 
trum  analysis.  'The  vertical  bar  indicates  the  velocity  which  the  manufacturer's  algo¬ 
rithm  calculated,  yet  there  are  two  peaks,  one  of  which  is  false  and  the  other  which  yields 
the  data  desired.  Unfortunately,  there  were  only  five  sets  of  spectra  with  this  feature 
which  were  obtained  during  the  first  month  of  operation,  therefore  only  a  subjective 
analysis  based  on  a  very  small  data  set  could  be  achieved.  By  comparing  the  velocity 
variance  contained  in  the  sLx-minute  data  to  five  sets  of  spectra,  a  velocity  variance  of 
3.0  m  5-'  was  chosen;  and  it  was  only  in  the  lowest  2  km  where  the  velocity  variance  ever 
exceeded  this  value,  possibly  due  to  extreme  turbulence  within  the  inversion  layer. 

If  each  of  these  criteria  are  applied  to  the  six-minute  profiles,  the  six-minute  profiles 
should  resemble  the  manufacturer's  consensus.  Fig.  31  shows  these  side-by-side  com¬ 
parisons  of  the  six-minute  data  and  the  manufacturer's  consensus.  Each  of  the  colors 
represent  data  which  would  be  thrown  out  by  the  three  thresholds  determined  earlier  in 
this  section:  red  barbs  are  winds  which  exceed  the  velocity  variance,  cyan  barbs  exceed 
the  power  differential  between  the  two  beams,  and  green  barbs  exceed  the  minimum 
power  requirements.  Black  wind  barbs  conversely  are  classified  as  the  good  data.  It 
appears  evident  from  this  final  product  that  quality  control  within  the  six-minute  data 
is  definitely  useful,  in  that  the  black-colored  barbs  match  well  with  the  consensus  data. 
Generally,  the  height  of  the  radar  is  limited  by  the  return  strength  of  the  signal,  w'hich 
is  related  to  the  amount  of  humidity  present.  Preliminary  data  from  Kadena  airfield  site 
in  Japan  during  the  Tropical  Cyclone  Motion  experiment  (TCM-90)  proves  that  higher 
heights  are  obtainable  in  a  deep,  moist  atmosphere,  as  seen  in  Fig.  32.  There  are  bad 
data  in  the  lower  atmosphere  as  well,  which  are  also  related  to  beam  difference  problems 
and  excessive  turbulence. 
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Plot  of  Various  Power  Difference  Thresholds  Between  Off-vertical 
Beams:  Combined  Heights  from  5197-9197  m  for  0605-1505  UTC  02 
May  1990. 
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Plot  of  Beam  Power  Difference  vs  Velocity  Standard  Deviation  and 
Mean:  Combined  Heights  below  3000  m  for  0605-1505  UTC  02  May 
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Fig.  30.  E.xaniple  Spectrum  Depicting  Velocity  Shift  due  to  Double  Peak 
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Profiles  from  02M8y9 


OaJo  from  J7-SEP  00:03  through  t7-SEP  23:01  UTC 
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Fig.  32.  Representative  Consensus  Data  from  TCM-90 
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C.  INVESTIGATION  OF  VARIABLES  AFFECTING  CONSENSUS 

The  two  variables  which  can  be  feasibly  changed  within  the  manufacturer's  consen¬ 
sus  procedure  (i.e.,  software  modifications)  are  the  minimum  number  of  points  accepted 
for  the  consensus  and  the  velocity  threshold.2  The  minimum  consensus  number  the¬ 
oretically  cannot  be  improved  because  the  occurrences  of  incomplete  hours  (satellite  in¬ 
terference  automatic  shutdown  periods,  power  outages,  frontal  passages,  etc.)  prohibits 
reducing  the  value  below  four.  This  is  dictated  by  research  done  by  NOAA's  Wave 
Propagation  Laboratory  (as  outlined  in  their  contract  specifications  to  purchase  a  net¬ 
work  of  wind  profilers  from  UNISYS).  Therefore,  an  attempt  was  made  to  improve 
consensus  by  varying  only  the  velocity  threshold  from  2.0  m  s-\ 

Subjectively,  Fig.  33  illustrates  that  by  varying  the  velocity  threshold  from  1.0  to 
5.0  msr'  there  generally  is  little  effect  on  the  final  hourly  consensus  profile.  However, 
when  the  velocity  threshold  is  diminished  below  1.0  m  s-\  there  are  noticeable  deviations 
regarding  direction  of  the  individual  wind  barbs  from  the  Tycho-provided  consensus, 
which  is  most  notable  below  2000  m.  To  validate  the  visual  conclusions,  the  difference 
between  the  new  consensus  (obtained  wth  the  new  threshold)  and  the  manufacturer's 
consensus  was  calculated.  Then  the  error,  or  "departure",  of  each  individual  six-minute 
datum  was  computed  by  subtracting  the  six-minute  velocity  from  the  new  consensus. 
Statistically,  differences  between  the  six-minute  u-velocity  and  v-velodity  components 
and  the  consensus  confirm  what  is  seen  graphically,  in  that  there  arc  numerous  "depar¬ 
tures"  greater  than  2.0  m  s-'.  Graphically,  there  arc  also  notable  deviations  from  the 
Tycho  hourly-consensus  when  a  velocity  threshold  of  5.0  m  s~'  is  exceeded,  which  are 
again  reflected  in  the  statistical  data.  For  a  velocity  threshold  above  5.0  m  s~'  there  are 
numerous  "departures"  of  the  sbe-minute  data  from  consensus,  i.e.,  "departures"  greater 
than  2.0  m  j"' 

These  results  indicate  that  there  remains  a  possibility  of  reducing  the  velocity 
threshold  below  2.0  m  s"',  in  order  to  improve  the  reliability  of  the  hourly  data.  This 
threshold  was  not  necessarily  important  for  the  period  on  which  this  study  concentrated, 
however,  in  other  atmospheric  conditions,  implementation  of  a  diminished  velocity 
threshold  may  be  useful.  For  example,  this  technique  of  reducing  the  velocity  threshold 
might  be  more  effective  higher  in  the  atmosphere,  where  the  signal  is  weak  and,  subse¬ 
quently,  the  consensus  number  would  be  low. 

2  Personal  conversation  with  Tycho  Industries'  technicians. 
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The  data  also  indicates  that  the  profiler  accuracy  can  fall  within  2.0  m  !r\  since  a 
consensus  can  be  obtained  with  the  velocity  threshold  set  at  1.0  ms“‘.  A  larger  data  set 
of  rawinsonde,  or  more  accurate  truth  data,  would  probably  confirm  the  accuracy  of  the 
profiler  to  be  within  1-2  m 
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Fig.  33.  Comparable  Hourly  Consensus  Profiles  from  02  May  1990.:  DilFcrcni 
hourly  consensus  profiles  are  plotted  utilizing  dificrcnt  velocity  thresholds 
varying  from  0.1  to  10.0  ms~K  The  hourly  consensus  based  on  the  man¬ 
ufacturer's  2.0  ntsr'  velocity  threshold  is  plotted  as  "Radar  Profiler  Con¬ 
sensus".  The  rawinsonde  data  associated  with  this  period  is  plotted  as  the 
"Sonde"  profile  in  the  left  column. 
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IV.  METEOROLOGICAL  ANALYSIS  USING  SIX-MINUTE  PROFILER 

DATA 


With  the  quality  of  six-minute  data  now  established,  as  well  as  some  limitations 
which  constrain  the  accuracy,  the  presence  of  high-frequency  meteorological  phenomena 
can  be  examined.  Reviewing  the  30  March  1990  to  05  May  1990  period  of  profiler  op¬ 
eration,  only  one  distinct  subperiod  in  the  pressure,  temperature,  wind  speed  and  direc¬ 
tion,  and  moisture  data  seemed  to  indicate  the  presence  of  any  high-frequency 
fluctuations.  The  subperiod  occurred  oh  02-03  May  1990.  Five-minute  meteorological 
data  from  03  May  1990  recorded  at  the  Monterey  Bay  Aquarium  Research  Institute 
(MBARI)  is  presented  in  Appendix  G.  The  wind  is  highly  variable  from  0000  to  0700 
UTC,  and  periodic  pressure  perturbations  occur  at  0800  and  1100  UTC  in  the  MBARI 
data.  The  pressure  perturbations  (approximately  1  to  1  1/2  mb  over  one-half  hour) 
though  subtle,  are  large  compared  to  the  normally  smooth  diurnal  change  experienced 
in  the  Monterey  Bay  area.  These  pressure  perturbations  have  been  verified  with  other 
high-temporal  resolution  barograph  data  from  the  NFS  Meteorological  Laboratory. 

A.  VELOCITY  AND  VELOCITY  PERTURBATION  FIELD  ANALYSIS 

* 

The  most  commom  display  of  six-nunute  profiler  data  is  the  "time  versus  height' 
field  of  velocity.  An  alternative  and  possibly  more  useful  way  to  examine  the  presence 
or  absence  of  high-frequency  phenomena  is  to  examine  the  velocity  perturbation  field; 
i.c,,  the  vector  velocities  remaining  after  a  mean  value  has  been  removed.  Both  of  these 
methods  will  be  utilized. 

Although  the  six-minute  data  of  horizontal  wind  components  can  yield  useful  and 
reliable  information,  detection  of  high-frequency,  low-amplitude  fluctuations  can  be  dif¬ 
ficult  due  to  the  relatively  high  mean  horizontal  wind  speed.  However,  vertical  motion 
measurements  could  reveal  high-frequency  perturbations,  since  the  mean  vertical  motion 
along  the  West  Coast  under  the  influence  of  the  subtropical  high  pressure  ridge  is  on  the 
order  of  centimeters  per  second.  In  fact,  periodic  motions  reveal  themselves  clearly  in 
Fig.  34,  a  time-height  series  of  vertical  motion.  A  dsmarcaticn  line  is  drawn  in  Fig.  34 
to  separate  general  regions  of  upward  and  downward  motions.  The  regions  are  charac¬ 
terized  by  maximum  vertical  velocities  on  the  order  of  1  m  s~K  The  demarcation  line 
appears  to  be  both  oscillatory  and  periodic,  approximately  2  hours  and  10  minutes. 
Because  the  approximate  structure  of  the  oscillations  are  too  small  to  be  sampled  by  si- 
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multaneous  measurements  from  the  available  weather  network,  some  principles  of 
time-space  conversion,  as  presented  by  Forbes  et  al.  (1989)  was  implemented.  Forbes 
states  that  "the  object  of  time-space  conversion  is  to  take  frequent  measurements  from 
a  single  location  while  the  prevailing  synoptic  winds  [and  high-frequency  propagating 
phenomena]  pass  overhead."  When  time-space  conversion  is  applied  to  Fig.  34,  an 
internal-gravity  or  Kelvin-Helmholtz-like  wave  is  suggested.  It  would  be  an  error  to 
assume  that  this  high-frequency  phenomena  is  southward  propagating  simply  because 
the  prevailing  winds  are  out  of  the  north.  Kelvin-Helmholtz  and  internal-gravity  waves 
traveling  along  a  stable  boundary  could  easily  propagate  against  the  prevailing  winds. 
More  information  from  other  sites  would  have  to  be  presented  to  determine  a  direction 
of  wave  propagation. 

Given  the  presence  of  some  periodic  forcing  in  the  vertical  motion,  it  is  anticipated 
that  corresponding  periodic  forcing  is  present  in  the  horizontal  velocities  as  well.  To 
examine  the  periodic  forcing  of  the  horizontal  velocity,  velocity  perturbations  were  cal¬ 
culated  by  subtracting  the  mean  velocity  from  the  individual  six-minute  observations 
(see  Fig.  35).  The  mean  was  calculated  at  each  level  from  the  long-term  consensus  for 
a  ten-hour  period  beginning  0600  UTC  02  May  1990.  The  horizontal  velocity  pertur¬ 
bation  components  reveal  a  similar  periodicity  as  noted  in  the  previous  paragraph  for 

the  vertical  velocity,  although  it  is  not  as  well  defined.  As  expected  from  the  "typical" 

* 

pattern  of  spring/carly-summer  winds  in  the  area,  there  was  the  presence  of  the  land  to 
sea  breeze  shift  (regions  highlighted  by  a  solid  line  at  approximately  2000  m  within 
Fig.  35).  However,  Fig.  35  shows  perturbed  winds  from  the  southwest  between 
4000-6000  m  from  0923  to  1012  UTC  and  from  the  northeast  from  1317  to  1406  UTC, 
which  correspond  to  two  of  the  periods  dominated  by  downward  motions  in  Fig.  34. 
These  regions  of  consistent,  large  perturbation  velocities  in  the  direction  of  the  prevailing 
winds  also  have  been  highlighted  in  Fig.  35  by  regions  bordered  by  thick  dashed  lines. 
Note  that  the  direction  of  perturbated  winds  in  these  latter  regions  are  opposite  in  di¬ 
rection  to  each  other,  yet  both  occur  during  corresponding  times  of  downward  motion 
in  Fig.  34,  which  is  hard  to  explain  given  the  limited  data  and  therefore  left  for  future 
study.  However,  the  directionality  of  these  perturbations  (aligned  northeast-southwest) 
suggest  that,  considering  conditions  over  the  ocean  are  uniform  and  the  perturbations 
are  due  only  to  the  high-frequency  source,  the  source  is  located  to  the  northeast,  or  in 
the  general  direction  of  the  pass  between  Atherton  Peak  and  Fremont  Peak. 
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Height  in  meters  (x  ld) 


Fig.  34.  Plot  of  NPS  Pronier  Vertical  Wind  Component  for  02  May  1990:  One 
inch  equals  1  m  j"'.  .Arrows  have  been  ofTsci  from  the  vertical  to  eliminate 
overlap. 
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.  35.  Plot  of  NPS  Profiler  Plorizontal  Wind  Perturbations  for  02  May  90:  One 
inch  equals  10  m  s~K 
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B.  THERMAL  ADVECTION  FROM  SIX-MINUTE  DATA 

To  better  understand  the  high-frequency  phenomena  observed  by  the  vertical  ve¬ 
locity  and  horizontal  perturbation  velocity  analysis  in  the  previous  section,  some  esti¬ 
mate  of  the  high-frequency  thermodynamic  properties  of  the  atmosphere  would  be 
useful.  Although  the  profiler  does  not  measure  thermal  properties  of  the  atmosphere 
without  a  special  acoustic  modification  (RASS),  some  aspects  of  the  thermal  structure 
can  be  inferred  from  profiler  wind  data.  In  particular,  a  single  wind  profiler  measures 
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From  these  quantities,  the  thermal  advcction  can  be  estimated,  given  some  rather  severe 
assumptions  about  the  (low. 

To  determine  the  thermal  advection,  it  is  assumed  that  all  Coriolis,  earth  curvature, 
frictional,  and  vertical  acceleration  terms  arc  negligible  compared  with  those  involving 
the  vertical  pressure  force  and  the  force  of  gravity.  This  leads  to  the  hydrostatic 
equation  presented  by  Holton,  1979 


where  clz  is  the  layer  thickness,  dp  the  pressure  difference  between  layers,  p  the  density, 
and  g  the  acceleration  of  gravity.  The  thermodynamic  state  of  the  atmosphere  at  any 
point  is  determined  by  the  values  of  pressure  (p),  virtual  temperature  (Tv,),  and  density 
(p)  at  that  point.  If  the  fluid  can  be  regarded  as  incompressible  these  field  variables  can 
be  related  to  each  other  by  the  equation  of  state  for  an  ideal  gas  as  presented  by 
Campbell  (1973) 

P  —  P^d^Vt  (2) 

where  Rj  is  the  gas  constant  for  dry  air.  Solving  Eq.  (2)  for  p  ,  and  then  substituting 
Eq.  (2)  into  Eq.  (1) 
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dz 


RdTy  dp 
g  P  ' 


(3) 


Integrating  Eq.  (3)  throughout  the  layer  gives 


fv.  =f 


dp 

g  P 


(4) 


The  term  on  the  left  side  of  Eq.  (4)  is  a  definite  integral,  and  Rj  and  g  are  constants; 
therefore,  from  Eq.  (4) 


(5) 


Next,  assume  that  temperature  is  linear  throughout  the  layer  and  Eq.  (5)  becomes 


(6) 


Two  more  equations  are  needed,  Holton  presents  the  hypsometric  equation,  which  re¬ 
lates  pressure  differences  to  height  or  thickness  changes  as 


Az  = 


g  Iv 


(7) 


Finally,  an  assumption  is  made  that  the  accuracy  of  the  geostrophic  wind  as  an  ap¬ 
proximation  to  the  observed  wind  is  a  good  one,  and  the  difference  is  negligible,  i.e., 
Coriolis  force  exactly  balances  the  pressure  gradient  force.  Gill  (1982)  presents  the 
geostrophic  approximation  in  component  form  as 

=  (8) 


g  dz 
f  dx' 


(9) 


or  in  the  vector  form, 
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(10) 


The  vector  difTerence  is  clearly 

aP,  =  (11) 

which  by  substitution  of  Eq.  (10)  into  Eq.  (I  I)  yields 

^  r  X  -  2,)  =  j  *  X  V;^Az).  (12) 

Substituting  Eq.  (7)  into  Eq.  (12) 

=  (13) 

but,  RJgl.  In(y7,//>j)]  =  A2/7',  by  Eq.  (7),  so  Eq.  (13)  now  becomes 

=  (14) 

By  definition  of  the  thermal  wind  by  Holton  (1972),  and  substitution  of  Eq,  (14)  into  the 
definition 


Finally,  temperature  advection  can  be  defined  as 


di 


TAd.  =  -y-^T,  or 
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(18) 


^  dT  dT 

^Adv  «  V  . 


Substitute  Eq.  (16)  and  Eq.  (17)  into  Eq.  (18),  and  replace  u  with  (the  consensus  ve¬ 
locity)  yields 


f^v 
Az  g 


]-v,C- 


fTy 
Az  S 


]. 


(19) 


Realistically,  and  are  not  known  so  u  and  v,  the  profiler-derived  winds,  are 
substituted.  Ty  was  calculated  from  a  formula  derived  from  Stull  (1988), 
Ty  =  {l  +  .6lq)T,  using  data  from  the  1300  UTC  03  May  1990  rawinsonde  launch. 
Large  deviations  c  way  from  geostrophic  balance  are  possible  for  six-minute  time  scales 
so  the  geostrophic  assumption  is  likely  to  produce  substantial  errors  in  V»VT.  Note 
that  this  analysis  will  highlight  regions  of  strong  vertical  shear  in  the  horizontal 
ageostrophic  (perturbation)  winds. 

When  temperature  advection  values  are  calculated  from  Eq.(19)  with  the  sbe-minute 
profiler  data,  the  result  is  generally  nothing  short  of  a  random  distribution  of  minimum 
and  maximum  values.  How’ever,  throughout  a  two-hour  period  on  03  May  1990,  de¬ 
picted  in  Fig.  36,  a  startling  consistent  pattern  of  concentrated  minimum  and  maximum 
values  lying  along  sloped  lines  is  revealed.  Again,  the  profiler  does  not  measure  tem¬ 
perature,  but  if  it  is  assumed  for  a  moment  that  Fig.  36  truly  is  temperature  advection, 
then  it  would  be  expected  to  also  have  packed  isotherms  in  the  regions  of  the  large 
advection  values  (i.e.,  packed  isotherms  would  be  needed  in  order  that  homogeneous 
atmospheric  vertical  motion  would  carry  "drastic"  changes  in  temperature.)  The  high 
values  of  temperature  advection  in  Fig.  36,  both  positive  and  negative,  are  probably 
associated  with  tighter  temperature  gradients,  but  not  directly.  The  slight  pressure  per¬ 
turbations  that  arc  induced  by  wave  activity  also  set-up  temperature  gradients,  similar 
to  the  packed  isotherms  Gossard  (1990)  revealed  regarding  turbulence  in  a  stratified  gas 
(sec  Fig.  37),  These  isotherm  gradients  in  turn  set  up  w'eak  velocity  perturbations  which 
are  subscqucnily  detccicd  by  ihe  profiler.  Although  this  phenomenon  is  easily  depicted 
using  this  temperature  advection  computation,  it  is  actually  only  an  artifact  created  by 
the  associated  weak  velocity  perturbations  observed  by  the  profiler. 

Subsequently,  the  vertical  velocity  for  the  period  represented  in  Fig.  36  was  also 
plotted,  as  Fig.  38.  Lines  indicating  regions  of  the  highest  density  of  extreme  temper¬ 
ature  advection  values  from  Fig.  36  were  transferred  to  Fig.  38.  There  is  a  tendency  for 


57 


(^01.  *]  sje^euj  U! 


Temperature  Advection  Field:  Reveals  areas  of  high  vertical  shear  In  the 
horizontal  winds. 


Fig.  37.  Packing  of  Isotherms  due  to  Turbulence  in  a  Stratified  Gas:  (Source: 

Gossard  (1990).) 

vertical  convergent  winds  (Fig.  38)  to  correspond  to  the  large  values  of  temperature 
advcction  (Fig.  36),  which  supports  the  idea  that  substantial  changes  in  the  horizontal 
velocity  perturbations  arc  occurring  across  these  layers. 

C.  SYNOPTIC  OVERVIEW  FOR  02-03  MAY  1990 

Synoptic  conditions  at  the  500  mb  and  surface,  representative  of  the  period  02-03 
May  1990,  are  presented  in  Fig.  39  and  40,  respectively.  The  combination  of  an  upper- 
level  low  height  center  over  southern  Arizona  and  a  high  pressure  ridge  along  the  West 
coast  force  prevailing  500  mb  winds  out  of  the  north  over  California.  Surface  conditions 
along  the  West  coast  arc  influenced  by  a  surface  pressure  trough  extending  from  the 
Arizona-New  Mexico  border  to  northern  California.  This  trough  creates  north  by 
north-westerly  flow  principally  parallel  to  the  California  coast;  however,  southerlies  ly¬ 
ing  along  the  California  central  coast  to  the  south  of  Monterey  are  against  the  synoptic 
scale  pressure  pattern.  Thi.s  suggests  a  possible  mesoscalc  eddy  or  low  pressure  center 
to  the  north  and  west  of  Monterey. 

Local  hourly-consensus  data  from  the  profiler  for  the  period  02-03  May  1990  are 
presented  in  Fig.  41  through  Fig.  44.  Winds  out  of  the  north  prevail  from  2000-6000 
m.  The  land-sea  breeze  shifts  are  evident  in  the  lowest  2000  m  of  the  atmosphere. 
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Vertical  Velocity  Field:  The  heavy  lines  represent  lines  of  large  temper¬ 
ature  advection  value  concentrations  transposed  from  Fig.  36. 


Specifically,  westerlies  depict  the  sea  breeze  crossing  the  Marina  shoreline  of  the 
Monterey  Bay,  which  persist  from  approximately  2100  UTC  to  0300  UTC.  From  0400 
UTC  to  2000  UTC,  the  land  breeze  sets  up  from  the  southeast,  very  much  aligned  to 
topography  of  the  Salinas  River  valley.  Hourly  observations  at  Salinas,  California  (20 
km  down  into  the  valley,  southeast  of  the  Monterey  Bay)  are  presented  in  Table  4.  The 
onset  of  the  land  breeze  is  evident  at  0500-0700  UTC  on  both  days  in  the  hourly  surface 
observations,  denoted  by  the  large  speed  decrease  in  the  northerlies. 

The  sbc-minute  data  for  this  period  of  land/sea  breeze  shift  is  presented  in  Fig.  45. 
As  seen  from  the  sbc-minute  data  presented  in  Fig.  45,  there  is  a  "burst"  of  easterlies  in 
the  lowest  4000  m  after  about  0414  UTC  02  May90  (as  well  as  on  03  May  90  at  ap¬ 
proximately  the  same  time).  (This  "burst"  is  more  evident  in  the  hourly  consensus  data. 
Fig.  41.)  The  easterlies  are  really  a  vector  summation  of  two  different  circulation  pat¬ 
terns:  (1)  the  land  breeze  from  the  south  by  southeast,  and  (2)  the  prevailing 
northeasterly  synoptic-scale  surface  flow.  This  "burst"  of  easterlies  occurs  because  the 
sea  breeze  (out  of  the  west)  dams  the  synoptic  northeasierlies  between  Atherton  and 
Fremont  Peaks,  and  when  the  sea  breeze  collapses  the  northcasterlies  now  flow  through 
the  pass  easily  and  abruptly.  (Refer  to  Fig.  2  for  topography  of  the  local  area.) 

D.  A  PROBABLE  EXPLANATION  OF  OBSERVED  HIGH-FREQUENCY 
PHENOMENA 

Although  the  determination  of  the  source  and  propagation  nature  of  the  oscillatory 
wave  seen  in  Fig.  34  arc  beyond  the  scope  of  this  study,  a  hypothesis  is  required  to  plan 
future  research  work.  Two  possible  explanations  are  that  some  type  of  internal-gravity 
or  Kclvin-Helmholtz  instability  waves  exist.  Although  the  wave  pattern  certainly  re¬ 
sembles  the  propagation  nature  of  an  internal-gravity  wave,  a  calculation  of  the 
Vaisala-Brunt  frequency  at  2000  to  4000  m  yields  a  value  or9.7  x  10-’ j-'-which  equates 
to  an  oscillatory  period  of  approximately  two  minutes.  This  period  differs  considerably 
from  the  two  hour  period  which  is  seen  Fig.  34.  Therefore,  1  conclude  that  the  high- 
frequency  phenomena  is  probably  due  to  shear  instability,  specifically  from  the  shear 
induced  by  the  reversal  of  the  sea  breeze  (i.e.,  land  breeze  southerlies  under  the  prevail¬ 
ing  northerlies  aloft).  This  is  particularly  convincing  as  the  oscillations  begin  shortly 
after  the  land  breeze  begins.  Additional  six-minute  profiler  data  can  be  used  to  confirm 
this  hypothesis  in  the  future. 
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Fig.  41.  NFS  Radar  ProFiIer  Hourly-consensus  Data  for  02  May  1990 
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Fig.  42.  NFS  Radar  Profiler  Hourly-consensus  Data  for  02  May 
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Fig.  43.  NFS  Radar  Profiler  Hourly-consensus  Data  fur  03  May  1990 
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Fig.  44.  NFS  Radar  Profiler  Hourly-consensus  Data  for  03  May  1990 
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Table  4.  SALINAS  HOURLY  SURFACE  WEATHER  OBSERVATIONS 
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V.  CONCLUSIONS 


This  study  determined  the  accuracy  of  the  NPS  profiler  to  be  3-4  ms-\  i.e.,  the 
standard  deviation  values  taken  from  the  atmospheric  layers  which  exceeded  a  .90 
corelation  coefficient  in  Chapter  III.  However,  by  obtaining  a  consensus  with  a  1  msr^ 
velocity  threshold,  which  was  similar  in  all  aspects  to  the  manufacturer's  consensus, 
suggests  that  the  accuracy  of  the  profiler  is  significantly  better  than  3-4  m  llius  there 
is  an  indication  that  the  profiler  does  meet  its  contractual  requirements.  This  disparity 
is  easily  explained  by  (I)  the  dislocation  of  sensors,  (2)  the  recent  evidence  ('Thomson 
and  Williams  1990)  that  rawinsondes  are  less  accurate  than  radar  profilers,  and  (3)  the 
extremely  small  data  set  available  for  this  study. 

Reliable  data  was  found  within  the  six-minute  profiles,  it  is  not  as  random  as  it  first 
appears.  The  parameter  thresholds  determined  in  Chapter  IV  are  very  effective  for 
identifying  bad  data  and  can  be  useful  in  improving  the  individual  sbt-minute  profiles. 
The  thresholds  also  can  be  used  individually  to  isolate  certain  physical  characteristics 
of  the  air  column  such  as  turbulence.  And  when  all  the  data  is  assimilated,  a  good  pic¬ 
ture  of  what  is  occurring  in  the  PBL  (up  to  6  km)  is  obtained.  Ironically,  even  the 
presence  of  bad  data  has  a  meaning,  such  as  the  presence  of  a  dry,  subsiding  air  mass. 
It  is  then  not  enough  to  say  the  data  is  bad,  it  must  be  understood  why  the  data  is  bad. 

Although  the  data  quality  checks  which  were  developed  to  clarify  the  six-minute 
data  were  very  successful,  the  presumption  that  they  wou^^  also  improve  the  hourly 
consensus  was  not  proven  with  as  much  success;  the  hourly  consensus  was  not  dramat¬ 
ically  affected  by  the  improvement  of  six-minute  data. 

Despite  the  large  vertical  spatial  resolution  (250  to  1000  m)  in  the  radar  data,  the 
overall  radar  performance  appears  to  be  very  good.  A  great  deal  of  detail  can  be  drawn 
from  the  wind  profiles.  It  has  been  showTi  that  high  frequency  meteorological  wave 
phenomena  can  be  detected  with  the  six-minute  data,  but  other  fluctuations  in  the  local 
meteorology  can  also  be  detected.  The  six-minute  data  can  accurately  depict  the  dra¬ 
matic  and  sudden  changes  in  the  land/sea  breeze  circulation.  Furthermore,  there  are 
complex  topographic  interactions  which  may  be  understood  by  analyzing  the  six-minute 
data.  The  high-temporal  resolution  and  accuracy  of  the  profiler  in  determining  three 
dimensional  meteorology  is  certainly  an  exciting  addition  to  the  NPS  Meteorology  cur¬ 
riculum. 
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VI. 


RECOMMENDATIONS 


This  study  has  been  able  to  demonstrate  the  quality  and  usefulness  of  the  profiler 
data  during  its  early  period  of  operation.  However,  there  is  considerable  potential  for 
future  and  more  complete  meteorological  studies.  The  following  list  of  recommen¬ 
dations  are  things  that  may  help  check  the  quality  of  the  data  and  better  define  the  de¬ 
tailed  meteorology  sampled  by  the  profiler. 

A.  MEASUREMENT  VALIDATION  AND  DATA  QUALITY  IMPROVEMENT 

•  Improve  the  manufacturer's  algorithm  for  spectral  averaging  based  on  local  mete¬ 
orology. 

•  Allot  a  rawinsonde  system  dedicated  to  work  with  the  profiler. 

•  As  Thomson  (1990)  indicates,  the  uncertainty  in  profiler-derived  \vTnd  speeds  and 
direction  is  clearly  significantly  less  than  the  natural  variability  of  the  •"'ind.  Rather 
than  using  rawinsondcs  which  are  highly  prone  to  sampling  errors,  another  radar 
profiler,  LIDAR,  or  SODAR  should  be  used  to  accurately  gage  NPS  profiler 
performance-not  the  rawinsonde  system. 

•  Investigate  techniques  such  as  those  listed  in  Appendix  D  which  can  improve  data 
quality  outside  the  hourly-profiles. 

B.  ACADEMIC  RESEARCH 

•  Use  a  50  MHz  profiler  for  comparison  with  the  NPS  profiler.  It  should  detect 
winds  higher  in  the  atmosphere. 

•  Investigate  the  possibility  of  acquiring  a  complete  network  of  at  least  three 
profilers,  either  through  outright  purchasing  or  coordination  with  other  agencies 
or  academic  institutions.  For  example,  San  Jose  State  University  might  be  inter- 
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ested  in  purchasing  a  profiler,  placing  it  within  the  Monterey  Bay  local,  and  then 
recover  the  data  remotely  via  telephone  lines. 

•  Add  the  NFS  2  kHz  SODAR  to  the  operational  data  base. 

•  Coordinate  5  kHz  SODAR  operations  with  the  Physics  Dept. 

•  Pursue  implementation  of  the  manufacturer's  software  regarding  Q,  calculations. 

•  Investigate  the  possibility  of  determining  a  Power  to  Humidity  algorithm. 

•  Rigorously  employ  the  Radio-Acoustic  Sounding  System  (RASS)  modification 
package  within  the  profiler. 

•  Make  all  remote  connections  to  automate  delivery  and  storage  of  the  profiler  data, 
especially  the  capability  to  upload  data  into  the  NPS  IDEA  Laboratory. 

•  Get  more  students  involved.  Draw  real-time  information  into  the  classroom. 

•  Conduct  a  detailed  study  of  land/sca  breeze  and  topography  interaction  along  the 
Monterey  Bay. 

•  Incorporate  real-time  ground-  and  space-based  remote  sensing  data,  such  as  wind 
profiles  from  the  NPS  Doppler  wind  profiler,  in  as  many  of  the  meteorology 
courses  as  is  possible. 
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APPENDIX  A.  EQUIPMENT  CONFIGURATION 
A.  WIND  PROFILER 


Fig.  A  -  I.  Functional  Diagram  of  Wind  Profiler 


'I'hc  wind  profiler  consists  of  a  two-stage  transmitter,  antenna,  receiver, 
transmit/receive  (T/R)  switch,  and  processor  with  its  firmware.  Fig  A- 1  illustrates  the 
basic  operation  of  a  Tycho-engineered  wind  profiler.  On  command  from  the  Real-Time 
Processor  (RTP)  the  antenna  controller  sets  up  the  proper  phasing  for  the  desired  an¬ 
tenna  beam  position  (i.e.,  vertical,  15®  southeast,  or  15®  southwest).  The  T/R  pulse  sets 
up  the  'f/R  switch  for  transmit  mode.  The  RTP  sends  a  prompt  signal  (or  pulse)  to  the 
receiver/modulator,  which  produces  a  pulsed  RF  signal.  This  is  then  disconnected  from 
the  antenna  during  the  transmit  cycle  to  prevent  overloading.  The  system  is  then  set  to 
the  "receive"  configuration.  The  receiver  is  reconnected  and  the  T/R  switch  directs  the 
weak  return  signal  from  the  antenna  to  the  receiver.  The  receiver  amplifies  the  signal 
and  extracts  the  in-phase  and  quadrature  (or  sine  and  cosine)  phase  components  from 
which  the  radial  velocity  is  derived.  Filtered  outputs  are  sent  to  the  RTP,  which  per¬ 
forms  the  various  processing  steps  necessary  to  produce  the  radial  wind  components, 
noise  levels,  signal  power  strength,  and  velocity  variance.  The  Profiler  Data  Handler 
(PDH)  is  a  VAX  3100  smart  terminal  which  processes  these  parameters  to  derive  time- 
height  wind  profiles,  Q  profiles,  wind  speed  and  direction  profiles,  radial  velocity  vari¬ 
ance  versus  height,  and  signal-to-noisc  (S/N)  profiles.  The  Keyboard  Display  Unit 
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(KDU)  has  four  functions:  (I)  it  provides  the  ability  to  edit  operating  parameters,  (2) 
displays  individual  mode  spectra,  (3)  monitors  profiler  operation,  and  (4)  controls 
profiler  operation. 

B.  RAWINSONDE 

The  upper  air  profiling  system  utilized  for  this  analysis  integrated  the 
Vaisala-DigiCORA  receiver-processor  with  Vaisala  RS80-series  rawinsondes.  The  spe¬ 
cific  software  version  installed  within  the  DigiCORA  system  was  the  Corex  version. 
Measurements  are  based  on  the  use  of  a  frce-ascent  balloon  rawinsonde.  The 
rawinsonde,  or  "sonde”,  transmits  data  to  the  station  at  402  MHz.  Refer  to  Fig.  A-2. 
Pressure,  temperature,  and  humidity  are  measured  by  sensors  in  the  sonde.  Each  of  the 
capacitive  sensors  controls  the  frequency  of  an  oscillator  through  an  electronic  switch. 
Relationship  between  the  frequency  and  corresponding  parameter  value  is  established  in 
the  calibration  process  prior  to  launch.  (Specifications  of  the  individual  sensors  are 
provided  in  Fig.  A-3.  The  location  of  the  sonde  is  determined  by  means  of  a  third  gen¬ 
eration  Navaid  (Omega  +  Very  Low  Frequency  (VLF))-based  windfinding  program 
called  "multifrequcncy  composite  VLF".  In  the  composite  solution  more  than  one  fre¬ 
quency  from  each  Omega  transmitter  and  one  frequency  of  the  USSR  VLF  network  are 
used.  Navaid  signals  are  relayed  to  the  ground  station  for  processing  and  for  wind  pro¬ 
file  computation.  Briefly,  the  signals  arc  processed  in  a  correlation  processor  to  derive 
their  relative  phases  from  which  the  sonde's  movements  can  be  calculated.  Wind  speed 
and  direction  over  the  range  of  0  to  180  ms-'  is  computed  with  an  accuracy  of  0.5 
ms-',  regardless  of  balloon  distance  from  the  ground  station. 


Fig.  A -2.  Vaisala  DigiCORA  Upper  Air  Sounding  System:  (Source:  Vaisala 


DigiCORA  Repair  Manual) 
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Table  A-1.  RAWINSONDE  METEOROLOGICAL  SENSOR  TECHNI 


CAL  DATA  (from  Vaisala  Corporation) 


Pnssure 

Type 

Capacitive  Aneroid 

Measuring  Range 

1060  hPa  to  3  hPa  (mb) 

Resolution 

0.1  hPa 

Accurao' 

0.5  hPa 

Temperaturt 

Sensor  Type 

Capacitive  Bead 

Measuring  Range 

+  60‘C  to  -90*C 

Resolution 

o.rc 

Accuracy 

0.2*C 

Humidity 

Sensor  Type 

HUMICAP*  Thin  Film  Capacitor 

Measuring  Range 

OtojOOVoR// 

Resolution 

\%RH 

Lag 

1  $  (6  m  J-‘  flow  at  1000  hPa,  +20“C) 

Accuracy 

2%  RH 

APPENDIX  B.  DERIVATION  OF  THE  RADAR  EQUATION  FOR  CLEAR 

AIR  PROFILERS 


The  radar  equation  relates  the  received  power  (echo  power)  to  radar  parameters  and 
target  parameters.  If  the  radiating  source  radiates  isotropically,  that  is,  equally  in  all 
directions,  the  power  density  (P, ,  power  per  unit  volume)  at  range  R  from  the  radiator 
is 


Pr  = 


(B-l) 


where  P,  is  the  peak  transmitted  power.  Radar  antennas  do  not  radiate  isotropically 
however,  but  rather  transmit  the  energy  in  a  beam.  The  energy  density  actually  received 
at  any  given  range  and  direction  relative  to  isotropic,  is  called  the  gain  (G)  or 

G  =  (P-2) 

Msotropic 

If  the  target  scatters  isotropically,  the  power  intercepted  by  the  receiver  is 

Pr  “  (iA  -X  ^ 


where  o  is  the  back-scatter  cross  section  of  the  target,  and  A,  is  the  effective  antenna 
aperature,  and  a  is  the  transmitter  antenna  cfilciency. 

The  widths  of  the  main  beam  in  the  north-south  and  east-west  directions,  6  and  <f> 
respectively,  are  generally  specified  by  the  w-idths  at  the  half-power  points.  These  widths 
depend  mostly  upon  the  size  of  the  antenna  and  wavelength  of  the  radiation,  although 
the  shape  of  the  antenna  and  power  distribution  across  the  antenna  also  influence  these 
widths.  To  a  good  approximation 

sin  d  =  ~ ,  {B  -  4a) 

sin  4>  =  ,  {B  -  Ab) 


Fig.  B-  1.  Illustratiun  of  Gconietrj'  for  Beam  Volume  Calculations 


where  L^  and  Lj  arc  the  linear  dimensions  of  the  antenna  in  the  directions  that  corre¬ 
spond  to  0  and  (f) ,  and  X  is  the  transnutted  wavelength.  For  narrow  beams,  the  sine  of 
the  angle  can  be  replaced  by  the  angle  itself: 

0==-~,  and  .  {B  -  5a) 


=  (B-Sb) 

■fhe  width  of  the  beam  between  the  first  nulls  (the  nulls  between  the  main  beam  and 
sidclobcs  on  either  side  of  the  main  beam)  is  about  twice  these  half-widths.  The  solid 
angle  corresponding  to  these  half-widths  (Q)  using  liqs.  (B-5a)  and  (B-5b)  is 

a  =  =  (/?-6) 


flic  area  of  a  rectangular  planar  array  (/!,)  by  substitution  into  Hq.  (B-h)  is  the  same 
as  the  physical  area: 


/f  =  L,/-2  = 


(/f-  7) 
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But  if  the  power  distribution  across  the  antenna  is  tapered,  the  effective  area  will  be 
slightly  different  from  this.  In  general,  the  effective  aperature  is  smaller  than  the  phys¬ 
ical  aperature,  the  ratio  being  called  the  antenna  efficiency  (a): 

a  =  ,or  Ag  =  ccA.  (B  —  8) 

The  antenna  gain  can  be  related  to  the  aperature  by  making  the  simplifying  approxi¬ 
mation  that  the  actual  beam  patiern  can  be  replaced  by  a  rectangular  beam  pattern  that 
has  a  given  gain  G  over  a  width  equal  to  the  half-width  and  zero  everywhere  else.  In  this 
case 


4nR^ 


but  by  substituting  Eq.  (B-7)  into  (B-9) 


Substituting  Eq.  (B-10)  into  Eq.  (B-3)  returned  power  becomes 


(B-9) 


(B-10) 


P,  P,G,A^fl,  .  (B  H) 

but,  a,  =  a,  =  a,  R,=  R,  (7,  =  (7,  =  (7,  and  A,  =  A,  =  A.,  because  the  antenna  is  the 
same  hardware  used  for  the  transmitter  and  receiver,  therefore 


TeTF- 

For  a  distributed  target  such  as  encountered  in  meteorology,  we  must  sum  over  all  the 
scatterers  in  the  sample  volume.  This  sample  volume  is  defined  by  the  radar  range  re¬ 
solution,  AR  =  ct/2  (where  x  is  the  pulse  length,  and  c  is  the  propagation  speed  of  radio 
energy),  and  by  the  half-widths  of  the  radar  beam  0  and  </>: 

2  2 
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Assume  the  scatters  don't  "shadow"  each  other  and  that  they  don't  have  appreciable 
absorption.  Then,  the  radar  reflectivity  is  the  sum  of  all  the  as  in  the  volume  divided 
by  the  volume: 


)7  =  -y-;  orcr  =  ?/K  (5-14) 

By  substitution  of  Eq.  (B-14)  into  Eq.  (B-12)  the  radar  equation  becomes 

Utilizing  Eq.  (B-10),  by  substitution  into  Eq.  (B-15) 

Utilizing  Eq.  (B-13),  Eq.  (B-16)  becomes 

P,  =  -2.6  «  p  -2~ .  (5-17) 

16rt^5^  87r5^ 

The  above  assumption  can  be  made  more  accurate  by  using  a  more  appropriate  form 
of  the  gain  equation.  If  a  Gaussian  beam  shape  is  used,  instead  of  a  rectangular  shape, 
the  fact  that  the  edges  of  the  sample  volume  arc  not  illuminated  as  powerfully  as  the 
central  portion  can  be  compensated  by  adding  a  shape  factor: 

5,  =  iP,A  Vex  — ~  ]  X  iShapeFactorl  (5-18) 

For  a  planar  array  according  to  Gossard  and  Strauch  (1983)  the  shape  factor  is 
4I9{X^IA).  Substituting  the  shape  factor  into  Eq.  (B-18)  the  radar  equation  evolves  into 

P,  =  P,AaV — (5-19) 
187r5^ 

Eq.  (B-19)  is  appropriate  for  the  scattered  signal.  The  radar  is,  of  course,  also  sensitive 
to  noise.  Noise  comes  from  a  variety  of  sources,  the  main  ones  being  system  noise  (also 
called  Johnson  noise)  and  cosmic  noise.  System  noise  can  be  controlled  somewhat  by 
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the  design  of  the  radar  but  it  can  never  be  completely  eliminated.  Noise  power  is  typi¬ 
cally  expressed  in  terms  of  the  equivalent  temperature  that  a  blackbody  would  emit  at 
the  given  frequency  to  account  for  the  power  observed.  The  noise  power  (?;,)  is 

P^=kBiaTc+Ts),  (5-20) 

where  k  is  Boltzmann's  constant,  B  is  the  receiver  bandwidth  (which  is  designed  to  match 
the  transmitted  bandwidth,  approximately  1/t  at  the  half-power  points),  Tc  is  the  cosmic 
noise  temperature,  and  is  the  system  noise  (also  called  the  receiver  noise).  The  factor 
a  arises  because  inefTiciencics  in  the  system  affect  the  noise  power  received  as  well  as  the 
signal  power: 


P,=  P,A,a'^cx^ 


n _ 

UR^k[pLTc+Ts)  ‘ 


(5-21) 


This  is  a  basic  equation  used  by  the  Tycho  profiler  and  is  consistent  with  research  done 
by  Rottger  and  Larsen  (1990).  Finally,  the  radar  constants  are  taken  out  of  the  equation 
and  grouped  together  for  ease,  as  used  in  the  body  of  the  study: 

p  ^  P,A,x^  Ar  ==  KP,a/  .  (5-22) 

Znk{aTc+Ts)  R^ 


where  is  a  "radar  constant". 
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APPENDIX  C.  RADAR  WIND  PROFILER  TRIGONOMETRY 


ZENITH 


Fig.  C-  1.  liorizont.'il  Velocity  Calculation  from  tlie  Three  Radial  Velocities 


Figure  C-l  illustrates  how  the  horizontal  wind  components  can  be  derived  from  the 
profiler-measured  radial  velocities.  Fig.  C-l  illustrates  three  I  'ams  which  arc  utilized  in 
either  the  low-  or  high-power  mode.  The  "zenith",  or  vertical  beam,  measures  the  w- 
component  of  the  wind  directly.  The  radial  velocities  of  the  other  two  ofT-vertical  beams, 
Xj^„  and  have  not  only  a  vertical  component,  but  aiso  horizontal  components 
related  by  the  following  equations: 

y  —  3'  '  4-  y  "  V  —  V  ’  ^  V  “ 

~~  ^RAD  ^  ^FjiO  ‘  RAD  ~  ‘  RAD  ^  *  RAD 

Xrad  ==  —  ( -  F;,  sin  0,  -f-  w  cos  Oj)  =  —  (w  cos  0,  +  u„  sin  0,) 


82 


Note:  The  negative  sign  compensates  for  the  scheme  depicting  velocities  toward  the 
radar  as  positive.  If  w,  and  (the  three  radial  velocities  processed  by  the 
profiler)  are  known;  then 

«/?  =  (-  YRad  -  vv  cos  0,)  /  shn  0,;  and 
v/j  =  -  ( -  -  w  cos  Oj)  /  sin 

Note:  A  negative  sign  must  be  applied  to  Vjj,  to  compensate  for  the  usual  horizontal 
cartesian  coordinate  system,  i.e.,  v  to  the  north.  0,  and  0^  are  the  angles  which  the  off- 
vertical  beams  are  offset  from  the  vertical.  These  angles  are  currently  set  at  15°.  ««  is 
the  Eastward  wind  component  computed  from  the  radial  velocities,  and  is  the 
Northward  wind  component  computed  form  the  profiler  radial  velocities. 
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w*/ -  upXcos(11.5*)  +  v,,x  sin(II.5‘) 

Component  1  Component  2 

Component  3  Component  4 

Fig.  C-2.  Trigonometry  Relating  Coordinate  Axes  in  the  Horizontal  Plane 
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Fig.  C-  3.  Trigonometry  of  Axis  Translation  for  Site  Correction 


site  Alignment 
Error  s  11.5® 
Site  Alignment 
I  55® 


Nots:  Correction  for  site  alignment  to  files  before  03  May  1990  2200  (Ul'C)  must  be 
applied.  After  03  May  1990  2200  (UTC),  the  elevation  and  orientation  arc  corrcct--no 
correction  needs  to  be  applied. 


84 


APPENDIX  D.  CONSENSUS  TECHNIQUES 

Because  we  know  that  spurious  data  can  enter  the  profiler  system,  much  more 
readily  that  a  direct-measuring  system,  automated  quality  control  (QC)  procedures  have 
been  developed.  This  appendbc  contains  excerpts  from  Brewster  (1989).  These  proce¬ 
dures  take  advantage  of  the  consistency  of  the  winds  in  time  and  height  (on  the  scale 
of  sample  separation,  250  meters  in  height  and  6  minutes  or  1  hour  in  time).  The  QC 
algorithms  for  the  profiler  data  calculate  an  expected  value  to  the  measured  value.  The 
expected  value  is  typicaliy  a  function  of  the  measurements  surrounding  the  datum  in 
time  and  space.  It  may  also  depend  on  some  sort  of  "first  guess"  provided  by  a  forecast 
model  or  expressed  as  a  function  of  the  surrounding  data.  If  the  difference  between  the 
expected  and  measured  values  exceeds  a  threshold,  then  the  datum  is  considered  suspect. 
Various  procedures  may  be  followed  to  determine  if  the  datum  arises  from  a  plausible 
but  extreme  meteorological  event  or  if  the  large  difference  is  caused  by  a  bad  datum 
among  those  used  to  calculate  the  expected  value.  In  this  appendix  several  techniques 
that  can  be  applied  to  wind  profiler  data  arc  described-first  that  which  is  developed  for 
use  by  the  Tycho  profiler,  then  some  others  that  can  be  applied  in-house  by  forecasters 
and  are  included  for  reference  as  techniques  to  apply  in  future  studies  using  the  NFS 
profiler. 

For  wind  profiler  data  distributed  to  data  networks  in  the  future  the  QC  must  be 
done  in  real  time,  that  is,  the  data  for  the  current  hour  must  be  checked  before  data  from 
the  next  hour  arrive.  Estimating  measured  values  is  more  difficult  because  the  data  are 
unevenly  distributed  in  time-plentiful  data  from  past  times,  but  none  from  future  times. 
For  some  applications,  more  precise  "two-sided"  quality  control,  done  after  future  data 
have  arrived,  is  not  only  possible  but  desirable.  Two  factors  merit  consideration:  1)  how 
long  after  the  data  collection  can  one  wait  before  running  the  QC,  and  2)  how  sensitive 
is  the  application  to  data  errors. 

A.  CONSENSUS  AVERAGING 

The  1-hour  wind  averages  produced  at  NFS  are  derived  from  the  6-minute  radial 
velocity  samples  by  a  technique  known  as  consensus  averaging.  Frimarily  an  averaging 
process,  this  technique  also  has  built-in  QC.  As  suggested  by  the  name,  not  all  the  data 
are  averaged,  but  only  those  that  mutually  agree  within  a  threshold.  A  2.0  m  s**  window 
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is  used  for  the  velocity.  We  review  the  effect  of  this  technique  by  considering  these  ten 
6-minute  radial  velocity  samples: 

3.3  4.6  20.9  2.8  3.6  4.1  -6.8  3.4  4.0  22.2  (ms“’) 

The  average  of  all  samples  is  6.2  m  s"'.  Looking  at  the  indi/idual  points  you  can  see  that 
6.2  m  s-'  does  not  at  all  represent  the  typical  sample  value  during  the  hour;  it  has  been 
severely  contaminated  by  three  outliers.  Now  if  we  collect  with  each  sample  all  other 
samples  agreeing  within  ±  1.5  m  s*' ,  that  is  falling  within  a  3.0  m  s-'  window,  we  get  the 
groups  from  Table  D-I.  The  greatest  number  of  samples  in  a  group  is  seven,  a  tie  among 
five  groups.  The  group  of  seven  coming  latest  in  the  hour  (indicated  by  the  arrow)  is 
then  selected  for  the  consensus;  in  this  case  all  five  groups  of  seven  are  identical  so  the 
position  within  the  hour  doesn't  matter.  The  consensus  average  is  the  average  of  that 
group,  3.7  m  s~\  which  almost  anyone  would  agree  is  representative  of  the  entire  sample 
set. 

Table  D-I.  FUNCTIONAL  CONSENSUS  PROCESS:  Example  of  individ- 
ual  sample  cross-comparisons  involved  with  consensus  process. 


sample 

samples  within  window 

3.3: 

2.8 

3.3 

3.4 

3.6 

4.0 

4.1 

4.6 

4.6: 

3.3 

3.4 

3.6 

4.0 

4.1 

4.6 

20.9: 

20.9 

22.2 

2.8: 

2.8 

3.3 

3.4 

3.6 

4.0 

4.1 

3.6: 

2.8 

3.3 

3.4 

3.6 

4.0 

4.1 

4.6 

4.1: 

2.8 

3.3 

3.4 

3.G 

4.0 

4.1 

4.6 

-6.8: 

-6.8 

3.4: 

2.8 

3.3 

3.4 

3.6 

4.0 

4,1 

4.6 

4.0: 

2.8 

3.3 

3.4 

3.6 

4.0 

4  1 

4.6  -< — 

22.2: 

20.9 

22.2 

In  practice,  mathematical  shortcuts  speed  calculation  of  the  average  value  that  has 
the  most  samples  lying  within  the  threshold  window.  If  no  group  has  at  least  four 
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members  to  form  a  consensus  (as  might  be  found  if  the  data  are  made  up  of  random 
'-’alues),  then  the  average  value  is  reported  as  zero. 

Consensus  averaging  handles  isolated  spikes  in  the  data  quite  well,  as  in  the  exam¬ 
ple.  It  also  handles  cases  in  which  the  returned  atmospheric  signal  is  too  weak  to  .be 
detected.  Then,  the  data  samples  are  nearly  random,  and  a  missing  value  will  be  re¬ 
ported. 

B.  COMBINED  MEDIAN  FILTER  AND  SHEAR  CHECK 

As  the  name  suggests,  this  routine  consists  of  two  checks  applied  sequentially  to 
hourly  averaged  data.  In  the  median  check  (see  Fig.  D-1)  data  are  gathered  from  adja¬ 
cent  hours  (the  current  hour  and  the  last  two  available  hours)  and  adjacent  levels.  The 
median  (the  reported  value  of  which  half  the  samples  are  greater  and  half  are  less)  of  this 
collection  of  data  is  then  computed  for  the  two  horizontal  wind  components.  If  the 
difference  between  the  datum  and  median  is  greater  than  a  threshold,  the  median  is  re¬ 
calculated  using  only  those  data  from  the  current  hour  and  the  last  hour.  This  is  done 
to  preserve  data  in  the  event  of  a  strong  wind  shift  with  time  as  might  be  found  in  a 
dramatic  trough  passage.  Data  are  flagged  bad  if  the  observed  datum  and  the  recalcu¬ 
lated  median  also  differ  by  more  than  the  threshold.  The  threshold  value  depends  on  the 
height  of  the  datum  (greater  thresholds  at  greater  altitudes),  the  wind  speed  (greater 
thresholds  for  high  speeds),  and  the  time  difference  between  the  datum  and  ihc  data 
collected  for  the  median  computation  (greater  threshold  when  there  arc  missing  data  in 
previous  hours).  In  the  current  implementation,  the  median  check  is  quite  loose;  it  is 
designed  to  flag  the  most  widely  varying  data  but  may  pass  along  a  substantial  number 
of  spurious  data.  The  philosophy  is  to  discard  the  implausible  values  before  the  shear 
check  and  to  leave  the  more  difficult  decisions  to  that  second  step. 

The  vertical  consistency  check  uses  a  gatc-to-gate  shear  threshold,  which  is  a  func¬ 
tion  of  height,  distance  between  the  observations,  wind  speed,  and  the  diflcrence  in  wind 
direction  (greater  .shear  is  allowed  if  the  wind  direction  is  unchanging).  A  typical 
threshold  is  about  20  ms-'  per  kilometer  or  about  13  kts  per  thousand  feet. 

If  the  magnitude  of  the  vector  shear  between  the  first  and  second  gates  exceeds  the 
threshold,  the  shear  between  the  second  and  third  gates  is  checked.  If  the  second  shear 
is  nearly  the  same  magnitude  as  the  first  and  in  the  opposite  direction,  the  data  in  the 
second  gate  are  flagged  bad  (both  horizontal  components).  This  situation  is  illustrated 
in  Fig.  D-2.  The  magnitude  of  the  shear  between  the  first  and  second  gates,  22.6  kt,  is 
excessive,  thus  prompting  a  check  on  the  shear  between  the  second  and  third  gates  (27.5 
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Fig.  D  —  1.  Data  Points  from  Mhich  Median  is  calculated  in  the  Median  Check 

ki).  The  two  shears  are  comparable  in  magnitude  and  nearly  opposite  in  direction  (the 
components  have  opposite  signs).  Thus  the  light  northwest  wind  in  the  second  gate  is 
flagged. 

If  the  wind  in  the  third  gate  tends  to  continue  the  shear  between  the  first  and  second 
gates,  a  straight  line  is  plotted  through  the  u-  and  v-components  of  winds  in  the  first  and 
fourth  gates.  If  both  components  in  the  second  and  third  gates  lie  reasonably  close  to 
their  respective  lines,  they  arc  accepted  as  valid  measurements.  This  is  the  condition 
shown  in  Fig.  D-2b,  where  the  sample  winds  imply  strong  shear  across  a  frontal 
boundary.  On  the  other  hand,  if  cither  component  in  the  second  and  third  gates  lies  far 
from  the  line,  both  winds  are  flagged  as  bad.  Fig.  D-2c  illustrates  the  case  of  two  spu¬ 
rious  winds  between  two  valid  winds.  The  line  for  the  u-components  lies  along  the 
horizontal  axis.  Although  the  middle  two  u-components  lie  close  to  the  line  for  u,  the 
middle  two  v-components  lie  far  from  the  line  for  v.  Thus  the  middle  two  winds  are 
flagged.  Note  that  the  scale  for  u  and  v  in  Figs.  D-2b  and  D-2c  has  been  held  constant 
to  facilitate  comparisons. 

'fhe  practical  result  of  these  tests  is  that  a  layer  of  winds  containing  four  or  five 
measurements  that  pass  the  median  check  and  corroborate  each  other  will  be  retained 
even  if  there  is  substantial  shear  at  the  layer  interfaces.  Moreover,  two  adjacent,  spuri¬ 
ous  winds  supporting  each  other  are  likely  to  be  flagged  correctly  as  bad,  as  in  the  above 
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Fig.  D-2.  Tiiree  examples  of  tlie  vertical  consistency  check.;  (a)  A  bad  wind 
bciween  two  good  winds,  (b)  good  winds  with  strong  shear,  and  (c)  bad 
winds  in  gates  2  and  3. 
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example;  three  such  winds  may  not  be.  Conversely,  if  a  layer  of  anomalous  but  valid 
winds  is  only  samples  by  two  or  three  gates,  those  wind  observations  may  be  erroneously 
rejected. 

The  vertical  consistency  check  is  applied  once  to  the  low  mode  data  and  then  to  the 
entire  profile,  high  and  low  modes  meshed  together.  Checking  in  each  case  begins  with 
the  lowest  gate  and  proceeds  upward.  Special  processing  with  linear  extrapolation  of 
neighboring  data  is  used  at  the  lowest  and  highest  gates  in  the  profile,  should  large  shear 
be  indicated  there. 

To  re-emphasize,  the  most  powerful  check  in  the  system  is  the  vertical  consistency 
check.  Substantial  changes  from  one  hour  to  the  next  are  generally  allowed  to  pass 
through  the  median  check  but  are  rejected  later  if  there  is  insufficient  continuity  in  the 
w'ind  shift  with  height. 

C.  OPTIMAL  INTERPOLATION 

Statistical  objective  analysis,  also  known  as  Optimal  Interpolation  (01),  is  widely 
used  for  analysis  of  meteorological  variables  on  quasi-horizontal  surfaces.  In  fact,  01 
is  used  today  at  many  national  weather  centers,  including  the  National  Meteorological 
Center  (NMC)  and  Fleet  Numerical  Oceanography  Center  (FNOC),  for  model  initial 
analysis.  This  method  utilizes  statistical  properties  of  the  observations  and  the  back¬ 
ground  field  (also  called  a  "first  guess").  The  background  is  usually  a  model  forecast  but 
climatological  averages  or  persistence  can  serve  the  same  purpose.  The  01  analysis  takes 
into  account  the  positions  of  the  observations  relative  to  each  other  and  relative  to  the 
analysis  point;  it  has  an  extrapolative  property  that  aids  analyses  near  data  edges. 

Optimal  interpolation  can  be  used  for  quality  control.  The  datum  under  scrutiny  is 
compared  with  a  value  analyzed  at  the  same  location  by  means  of  01.  Surrounding 
observations  (but  not  the  central  one)  and  the  background  field  contribute  to  the  ana¬ 
lyzed  value.  If  the  difference  between  the  analyzed  and  observed  value  exceeds  some 
threshold,  the  datum  requires  further  testing.  Since  the  01  method  provides  an  estimate 
of  the  analysis  error,  the  threshold  can  be  specified  as  a  function  of  the  analysis  error. 
This  allows  for  more  variation  when  observations  arc  sparse  or  distant  from  the  datum 
in  question  and  less  variation  when  the  observations  arc  dense  or  located  very  close  to 
the  datum. 

Large  discrepancies  between  the  observed  and  analj’zed  value  may  arise  either  be¬ 
cause  the  subject  obser%'ation  is  erroneous  or  because  one  of  the  observations  contrib¬ 
uting  to  the  analysis  is  erroneous.  To  find  out  which  is  true,  the  data  are  re-analyzed 
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several  times,  removing  one  of  the  surrounding  observations  each  time  to  see  if  that 
produces  agreement  (Fig.  3).  The  removed  observation  is  returned  for  subsequent  re¬ 
analyses  so  that  only  one  observation  is  absent  for  each  re-analysis.  If  no  agreement  is 
found  for  any  re-analysis,  the  datum  being  examined  is  flagged.  Should  the  removal- of 
one  of  the  observations  produce  agreement,  the  removed  observation  is  flagged  as  sus¬ 
pect  and  is  not  used  in  any  subsequent  QC  analysis.  Suspect  flags  are  internal  to  the 
QC  program  and  are  generally  not  part  of  the  QC  output.  Typically,  a  datum  is  flagged 
only  when  it  is  the  subject  datum. 

This  technique  just  described  can  be  used  to  check  data  in  a  time-height  cross- 
section  as  well  as  on  the  horizontal  surfaces  for  which  it  was  originally  designed.  The 
correlation  statistics  must  then  be  specified  in  time-height  space.  01  time-height  quality 
control  of  profiler  data  in  real  time  presents  a  special  problem.  Because  there  Is  an  at¬ 
tempt  to  always  examine  the  latest  hour,  the  analyses  are  not  centered  in  time  (i.e.,  it  is 
always  analyzing  on  the  edge  of  the  time-height  domain).  A  second,  often  vexing 
problem  is  that  spurious  winds  are  not  always  isolated-repeated  errors  in  a  given  beam, 
extending  over  several  gates,  can  cause  the  rc-analysis  technique  to  fail. 

The  majority  of  cases  can  be  handled  properly  if  data  selection  for  each  analysis  is 
done  to  minimize  the  possibility  of  multiple  bad  datj,  points  being  included.  An  01  QC 
algorithm  developed  at  PROFS  makes  two  passes  through  the  data  with  the  flagging 
threshold  decreasing  on  the  second  pass.  In  the  first  pass,  data  are  checked  in  sequence 
from  the  lowest  gate  to  the  highest.  The  two  nearest  data  points  arc  collected  from  each 
of  the  numbered  sectors  shown  in  Fig.  D-4,  except  those  which  have  not  yet  been 
checked  (those  in  the  shaded  area-Scctor  I  and  pan  of  the  sector  containing  the  subject 
point  P).  The  eight  data  points  which  correlate  best  with  the  point  being  exanained  (as 
determined  from  historical  data)  are  used  in  the  analysis.  The  flagging  system  is  the 
same  as  described  in  the  paragraph  above.  The  second  pass  proceeds  as  the  first,  but 
data  are  accepted  from  all  the  sectors.  In  this  pass  it  is  presumed  that  if  there  are  any 
spurious  data  above  the  point  being  examined  (i.c.,  not  yet  checked  on  this  pass),  they 
have  been  flagged  bad  in  the  first  pass  or  flagged  suspect  when  they  w'erc  used  as  a 
contributing  observation  on  the  second  pass. 

D.  RECURSIVE  FILTER 

Quality  control  can  be  built  into  the  analysis  itself  by  adjusting  the  observation 
weights  based  on  the  observation  diflTcrence  from  a  first  quess  or,  more  commonly,  from 
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Fig.  D- 
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Continue  eliminating  each 
successive  observation 
while  retaining  all  the  others 


3.  Sclienialic  of  Successive  Eiiinitialiou  5’roce.ss 


92 


Fig.  D-'l.  Data  Search  Strategy  for  Optimal  Interpolation  Analysis  (01) 
QC:  Sector  searching  assures  an  even  distribution  of  data  around  the  da¬ 
tum  being  checked.  Data  within  shaded  area  arc  excluded  in  the  first  pass; 
all  sectors  arc  searched  in  the  second  pass. 

the  presious  analysis  in  a  multiple  pass  (successive  correction)  anal) sis  method.  In  this 
approach,  data  points  arc  not  discarded  but  arc  given  a  very  small  weight. 

The  recursive  filter  analysis  uses  bilinear  interpolation  to  place  the  obscr\cd  data 
(actuallv  the  difrcrcncc  from  a  first  guess  field)  onto  a  rcgularl)  spaced  grid.  To  smooth 
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the  analysis,  a  numerical  filter  (similar  to  a  weighted  running  average)  is  applied  to  the 
grid-point  values.  Then  the  process  is  repeated  using  the  last  analysis  as  a  first  guess. 
With  each  pass,  the  filter  is  modified  so  that  less  smoothing  is  done.  Quality  control  is 
introduced  by  weighting  each  observation  according  to  its  difierence  from  the  previous 
analysis.  One  can  discard  an  observation  if  its  "quality  weight"  becomes  very  small,  or, 
if  the  data  density  is  high  enough,  one  can  presume  that  observations  with  very  small 
weight  wDl  be  overwhelmed  by  those  with  large  weight  (good  quality  observations). 
Because  the  recursive  filter  requires  far  less  computation  than  optimal  interpolation  and 
no  data  searching  it  is  much  faster  than  01.  The  speed  depends  upon  the  number  of 
re-analysis  passes  necessary  and  the  number  of  filter  passes  performed  for  each  analysis 
pass.  A  good  first  guess  can  make  the  number  of  analysis  passes  very  small  (e.g.,  five) 
and  can  greatly  increase  the  accuracy  of  the  quality  control  decisions. 

This  technique  has  been  tested  on  six-minute  profiler  data  at  PROFS.  The  analysis 
is  done  on  a  lOO-m  by  6-minutc  grid  and  the  first  guess  consists  of  the  median  of  11 
observations  around  each  grid  point.  A  time-to-space  conversion  which  equates  six 
minutes  to  100  meters  is  used,  and  data  are  discarded  (QC  flag  set  to  "bad")  if  their 
weight  is  less  that  lO'h 
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APPENDIX  E.  INDIVIDUAL  PROFILER  COMPARISONS  AND 
RAWINSONDE  LAUNCHINGS 
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Fig.  E-  I.  Profiler-RaHinsonde  Wind  Coinparisoiis  for  20  April  1990 
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Fig.  E-2.  R.'iu  insonde  Soimding  from  20  April  1990 
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E“3.  Profiler- RaMinsonde  Wind  Comparisons  for  1600  UTC  23  April  1990 
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Launch  time  ;  90  04  23 
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E-4.  Ra>vinson(le  Sounding  from  1600  UTC  23  April  1990 
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Fig.  E-7.  Profiler- Raw  insonde  Wind  Comparisons  for  2(K)0  UTC  02  May  1990 
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Launch  time  :  90  05  02 
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E-8.  Ra\un.soiide  Sounding  from  1900  UTC  02  May  1990 
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Fig.  E-9.  Prolller-Rauiiisomle  Wind  Comparisons  for  0100  UTC  03  May  1990 
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E-  to.  RaMinsonde  Sounding  from  0100  UTC  03  May  1990 
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Temperature  {"C)  Speed  (m/s)  Direction 


Launch  time:  90  05  05  at  6  17  GMT.  Sil.e,  Lot:  36.300,  Long:  -122.48 
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Pig.  E  -  II,  Profiler- Rauinsoiule  Wind  Comparisons  for  0700  UTC  03  May  I9‘I0 
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Fig.  E-  12.  Ra^unsontle  Sounding  from  0700  UTC  03  May  1990 
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Fig.  E-  13.  Profilcr-RaMiiisomle  Wind  Coinp.nrisoijs  for  13(10  UTC  03  M.iy  1990 
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Sife  Lot:  36.300.  Long;  -122.800 
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E-  [‘j.  RaHiiisonde  .Sounding  from  13UI)  UTC  03  May  1990 
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Profiler  —  Rawinsonde  Comparison  Winds 

Launch  iime:  90  05  05  at  18  26  GMT.  Site  Lar:  36.250,  Long:  123.23 
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Fig.  E-  15.  Profiler- Ra« insonde  Wind  Coinp.msuns  for  19U0  UTC  03  May  1990 
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Fig.  E-  16.  Rawla.sonde  Sounding  Trom  1900  UTC  03  May  1990 
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APPENDIX  F,  NFS  PROFILER  SITE  SURVEY  RECORD  AND  LAYOUT 


Naval  Poot{)i'«duate  School 
Depat'ttiieril:  of  Oceariof)i'aphy 

i 

*  19  OPfl  90 


Frofiii  CoiiUitancJot*  Kui't  Schncbslii,  NOAD  (Codo  ClC/Sn) 

Toi  Proruuuor  UilliaM  Shaw  (Code  MIU 

Sub)i  nZIHUlM  t  POSITION  OF  FORT  ORD  PROFILER 

End  I  (1)  nzitiuith  t  P'osition  Dotartiiinat  ion  of  Foi't  Drd  Profile)* 
<li>  Fo»*t  Oi'd  Pi'ofllBi'  Survey  Uorhaheuta 
<31  Pouitions  t  Oeacrlptions  of  Control  Stations 


I.  At  the  r{?(|ueBt  of  LCOR  Oan  Strucd,  USN;  Hi*.  Jariins  Cherry  and 
1  completed  a  survey  to  deterMine  a  peodetic  aziiiiuth  arid  poultiori 
foi*  tlio  pi'ofilor  array  at  Fo»*t  Ord.  Thu  results  are  aumiiiai'i  sod 
below  and  dlocuuned  in  Enclosure  (I). 


OslMiith  of  Side  336*  07.0*  (linear  error  SSX  a  3.3*  > 
(wuuturn  side)  • 

Asihuitti  of  North  021*  07.0’  (auauiiiing  array  Is  square) 
South  Dlanonal 


0i*ray  Center 
Position  on 
NAD  03  datuhi 


N  3G»  41’  A3. 716"  (circular  ori'or  95v!) 
W  121*  AS*  29. oaA"  (a  3  cm  relative  to) 
(station  DAIRY.  I 


Elevation  SI  Miiters  above  sea  level  (NOVO  29) 

Ai*ray  Top 


2.  Note  that  the  elevation  was  duteriiiined  froM  only  one  point 
(station  DAIRY).  Lachinp  an  independent  cliech  fi’OM  anothui* 
survey  point,  it  is  best  to  auuuMo  that  elevation  is  only  n<^'C'd  to 
±1  Mutur.  Doth  the  position  and  azlMuth  were  dutei'Mlned  ustn,'| 
r'udundant  MeauureMunts  in  a  1 eaut-squareu  adJustMunt,  whicli  nave 
tlio  precision  eutiMatuu  cited  above. 

3.  Note  that  we  sot  a  new  survey  MOnu)>ient  (DAIRY  OFFSET)  on  the 
hill  ovo.'loohltia  the  profiler.  If  needed,  the  axiMiith  could  bo 
clioched  froM  there  in  the  future. 


Kurt  J.  Schnebels 
CoMiiiandoi',  NOAA 
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MARIN  8  to  AIRPORT  BEACON 

rRlTZCIIE  AAP  PT  URD 


Fig.  F-1.  Atmospheric  Profiler,  Fort  Ord  Survey  Observ.ntiou  Sketch:  ('on- 
ductcd  19  April  1990  by  Cdr.  Schncbclc  (no  scale). 
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(deg  C)  (degrees) 


yVPPENDIX  G.  HIGH-FREQUENCY  METEOROLOGICAL  DATA 


Fig.  C-I.  Recorded  Data  for  03  May  1990.:  From  the  Monterey  Hay  Aquarium 
Research  Institute. 
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